Free energy of colloidal particles at the surface of sessile drops 
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The influence of finite system size on the free energy of a spherical particle floating at the surface 
of a sessile droplet is studied both analytically and numerically. In the special case that the contact 
angle at the substrate equals 7r/2 a capillary analogue of the method of images is applied in order 
to calculate small deformations of the droplet shape if an external force is applied to the particle. 
The type of boundary conditions for the droplet shape at the substrate determines the sign of the 
capillary monopole associated with the image particle. Therefore, the free energy of the particle, 
which is proportional to the interaction energy of the original particle with its image, can be of 
either sign, too. The analytic solutions, given by the Green's function of the capillary equation, 
are constructed such that the condition of the forces acting on the droplet being balanced and of 
the volume constraint are fulfilled. Besides the known phenomena of attraction of a particle to a 
free contact line and repulsion from a pinned one, we observe a local free energy minimum for the 
particle being located at the drop apex or at an intermediate angle, respectively. This peculiarity 
can be traced back to a non-monotonic behavior of the Green's function, which reflects the interplay 
between the deformations of the droplet shape and the volume constraint. 

PACS numbers: 68.03.Cd, 47.85.-g, 89.90.+n, 83.80.Hj 



I. INTRODUCTION. 

Two-dimensional (2D) structures formed by micron- 
sized colloidal particles trapped at fluid interfaces have 
received a significant attention since an interfacial col- 
loidal crystal has been first observed by Pieranski in 
1980 PJ. In the last two decades 2D-colloidal structures 
have proven to be an excellent model system for study- 
ing the influences of dimensionality and confinement on 
phase transitions. For example, a general scenario for 2D- 
melting, predicted by Kosterlitz and Thouless Q in 1973, 
has been directly verified in experiments using quasi-2D 
colloidal structures Q . It has also been found that strong 
lateral confinement can lead to peculiarities like reentrant 
freezing [J] . From the theoretical point of view the capil- 
lary interactions between colloidal particles mediated by 
a fluid interface resemble 2D-electrostatics, which is due 
to the fact that the small deformations of the interface 
are governed by the Poisson equation with the pressure 
field playing the role of a two-dimensional charge den- 
sity. In terms of application perspectives, controlling 
self-assembly [5| and structure formation @-[g at fluid 
interfaces provides tools for the engineering of advanced 
materials based on transferring such microstructures onto 
solid substrates |7J . Colloidal particles are also known to 
stabilize emulsions by self-assembling at the droplet sur- 
faces [9j, creatin g c apsules which might have potential 
use in medicine [10j . The particles at droplet surfaces 
can also be used to study correlations [IJ [l2[ , topolog- 
ical defects, and the generalized Thomson problem for 
optimal packing on a sphere [l3J. 

The presence of an interface strongly modifies the ef- 



fective interactions between colloidal particles as com- 
pared with their counterparts in the bulk. It has been 
already noted by Pieranski [l[, that in the generic case 
of different dielectric constants of the adjacent fluids, 
the interface breaks the symmetry in the ion distribu- 
tion around each particle, leading to an effective dipolar 
repulsion. Therefore, in the presence of lateral confine- 
ment a 2D-crystal can form. However, in the last decade 
many authors reported on the spontaneous formation of 
mesostructures like clusters, stripes (fJ-Qjj]: and freely 
floating crystallites [18| without confinement. These ob- 
servations suggest more complicated forms of the effec- 
tive interaction potentials, in particular the presence of 
an attraction, the origin of which has recently been in- 
tensively discussed in the literature. The interface medi- 
ated capillary interactions arc considered to be promising 
candidates for an explanation, but their theoretical de- 
scription is still in progress. The main idea, formulated 
already in 1949 by Nicolson [19J], is that the surface en- 
ergy of a deformed interface depends on the separation 
of the particles generating the deformation [20|. Origi- 
nally, Nicolson studied the behavior of millimeter-sized 
bubbles deforming the interface due to their buoyancy. 
In the case of micrometer-sized or smaller particles grav- 
ity can be usually neglected, but even then capillary in- 
teractions can dominate thermal energies by orders of 
magnitude [2l|, [22| ■ There seem to be at least two rea- 
sons for the origin of the interface deformation around 
the particles. First, in the case of charged colloidal par- 
ticles pinned at the interface, an electrostatic pressure 
emerges due to the inhomogeneous ion distribution close 
to the interface IH, 01. Second, the interface must de- 



form in order to meet the particle surface at a given con- 
tact angle or at a prescribed (pinned) contact line, which 
becomes important in the cases of smooth non-spherical 
particles [22], (25( or rough particles [21[ , respectively. 

The case of charged particles has been strongly de- 
bated in the last few years and different theoretical ap- 
proaches have lead to contradictory results for the form 
of the capillary interaction potential. In 2002 Nikolaides 
et al. |18[ observed a freely floating hexagonal structure 
formed by charged particles at the surface of a water 
droplet in oil. The authors proposed a simple explana- 
tion according to which the particles are pulled towards 
the water phase by the electrostatic forces due to the 
different dielectric constants of water and oil. They con- 
cluded that the presence of such a pulling force leads to 
a logarithmically varying interface profile around a sin- 
gle particle and consequently to a logarithmically vary- 
ing effective interaction potential between two particles. 
However, it was noted by Megens et al. [26| that this 
long-ranged interface deformation yields an unbalanced 
force acting on the rim of the vessel containing the liquid, 
which is inconsistent with Newton's third law. Instead, 
in the presence of the electrostatic pressure, decaying as 
r~ 6 with lateral distance r from the particle, the latter 
authors concluded that the interaction potential should 
be attractive and should decay also like r~ 6 . A similar 
reasoning was presented by Foret and Wiirger [23|, who 
obtained the same power law, but with opposite sign. 
Their calculation was based on a regularization of the in- 
terface deformation field at r — and a superposition ap- 
proximation for the electrostatic pressure in the presence 
of two particles. However, Domingucz et al. noted [241 ] 
that the latter assumption is not correct because the elec- 
trostatic pressure, proportional to the electrostatic field 
squared, is not additive. Instead, one should use the ad- 
ditivity of the electrostatic potential, which finally gives 
an attractive r -3 behavior. This leads to the conclu- 
sion, that the total interparticlc potential including the 
direct dipolar electrostatic repulsion, governed also by 
the power law r -3 , cannot exhibit a minimum, unless 
the screening length in water is comparable with the size 
of the particles. In such a case a shallow minimum with 
a depth of several ksT could be expected [2J|. However, 
its occurrence would depend sensitively on the precise 
values of several system parameters. 

In yet another publication [281 ]. Dominguez et al. sug- 
gested that the finite size of the system might be impor- 
tant for generating the long-ranged interactions. In fact, 
in the experiment reported by Nikolaides et al., the par- 
ticles were trapped at the surface of a droplet pinned to 
a solid plate, but curvature and volume constraint were 
neglected in the theoretical description in Rcf. [18| . For 
recent assessments of these capillary forces see Rcfs. [29J ]. 
[30| . and [31j . In the present work we show that finite 
size effects indeed lead to effects not present in the case of 
a flat unbounded interface. Quite generally, the presence 
of boundaries in a colloidal system can change the nature 
of the effective interactions. For example, the hydrody- 



namic interactions in the bulk decay like the inverse of 
the particle separation, but in the presence of one or two 
walls or of surfaces the decay exponent can change or 
there might be even a switch to a logarithmic law |32j. 
On the other hand, in one-dimensional channels the hy- 
drodynamic interactions are screened and characterized 
by an exponential decay [33|- Our ultimate objective is 
to study how curvature and confinement of the fluid in- 
terface affect the capillary interactions between the float- 
ing particles. The effects of curvature have been already 
studied by Wiirger [3J, ]35| , but this description cannot 
be regarded as complete as long as the consequences of 
force balance on the droplet have not been correctly re- 
solved [36| . In Ref . [37| Kralchevsky et al. have provided 
numerical results for the free energy of particles protrud- 
ing from a spherical liquid film and in some cases they 
observed a non-monotonic dependence of the ensuing ef- 
fective interaction on the spatial distance between them, 
but the physical mechanism responsible for this effect has 
not been explained. Oettel et al. [23| derived the equa- 
tion for the axisymmetric shape of a slightly deformed 
spherical droplet, taking into account the volume con- 
straint and a certain type of boundary condition at the 
substrate. However, the influence of the type of bound- 
ary condition on the free energy and the occurrence of 
a possible confining capillary potential due to the pres- 
ence of the substrate have not yet been discussed. In the 
following we aim at calculating such an effective poten- 
tial and at providing some general instructions of how 
the finite size effects should be taken into account when 
discussing capillary interactions. 

As a first step in understanding the interactions, we 
have performed a detailed analysis of the one-body prob- 
lem of a spherical particle at the surface of a sessile 
droplet (Sec. II). In Sects. Ill and IV we propose a 
perturbation theory for a slightly deformed droplet and 
solve the boundary value problem at the substrate by the 
method of images known from electrostatics. In Sec. V 
we present the results of the full numerical minimization 
of the free energy and compare them with the analytical 
results. Together with our conclusions we summarize our 
findings in Sec. VI. Various important analytic calcula- 
tions are presented in Appendices A, B, and C. 



II. MODEL 



We consider a smooth, solid spherical particle of radius 
a trapped at the surface of a liquid droplet of volume Vi 
and surface tension 7 residing on a planar substrate. For 
convenience we call the surrounding medium " gas" , but 
we assume a non- volatile liquid. The equilibrium contact 
angles 9q at the substrate and 9 P at the particle are given 



by Young's law: 
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where ^) a b are surface tensions with the indices a, b = 
0,p,l,g standing for substrate, particle, liquid, and gas, 
respectively, and "fi g = 7. In polar coordinates, the po- 
sition of the particle can be described by the radial dis- 
placement h and the polar angle a (see Fig. Q]). 

It can be shown [38| that in the absence of external 
forces the particle adjusts its immersion, independently 
of its lateral position at the droplet, according to Young's 
law (Eq. @) such that the droplet remains an unde- 
formed spherical cap. This configuration with the parti- 
cle at the unperturbed droplet will be called the reference 
configuration, for which we set h = 0. As a consequence 
of the spherical shape of the droplet the corresponding 
free energy does not depend on the polar angle a (see 
Fig. [!}. The only exceptions are configurations with a 
particle close to the contact line, where the movement of 
the particle is constrained by the substrate. This situa- 
tion is common in experiments with evaporating droplets 
due to the flux of liquid dragging the particles towards 
the contact line [331. When the height of the interface 
in the liquid wedge reduces below the particle diame- 
ter the interface around the particle must deform. This 
particular case could be discussed within the framework 
presented here, too, but it is beyond the scope of the 
present study. 



A. Free energy functional 

If, however, an external force / acts on the particle 
in radial direction, the interface deforms such that the 
capillary force counterbalances the external force. The 
corresponding equilibrium shape of the droplet minimizes 
the following free energy functional: 

•F[{r(fi)}, h, a; /, { 7a& }, a, V h A] = 

= l{Slg — Slg^ re f) + ("fQl — "fOg)(Sol ~ Sol,ref) 
+ (ipl ~ -fpg)(Spl ~ S pLre f) -fh- \{V - V t ), (3) 

where {r(il)} denotes the configuration of the interface 
in terms of spherical coordinates ft = (6, (p) on the unit 
sphere, and Si g , Sqi, S p i are the liquid-gas, substrate- 
liquid, and particle-liquid surface areas, respectively. In 
Eq. ([3]) T is the free energy functional relative to the ref- 
erence configuration of the unperturbed droplet referred 
to by the index ref. The last term with the Lagrange 
multiplier A ensures that the liquid volume is constant. 
The sign of this term has been chosen such that the value 
of A minimizing the functional in Eq. ([3]) can be inter- 
preted as the internal pressure. The sizes of the droplet 




FIG. 1: Sketch of the system (see main text); (a) refer- 
ence configuration; (6) deformation due to an external force 
/ = fe r for a free contact line at the substrate, which is 
not necessarily a circle but it is characterized by the an- 
gle d c (4>) as a function of the azimuthal angle <f> in the xy 
plane. Neglecting the line tension, in this case the contact 
angle at the substrate equals 6q along the contact line. The 
direction of the x-axis is chosen such that the center of the 
particle at ri, re / + he r (with h — in (a)) lies in the xz 
plane. Accordingly, the external force acts in the direction 
Qi — (0 — a, cj> — 0). In the case of a pinned contact line at 
the substrate (not shown), for / 7^ and a/0 the contact 
angle at the substrate varies along the contact line (i.e., it is 
a function of <j>) even if the contact line forms a circle. The 
contact line at the surface of the colloidal particle is taken to 
be free such that the corresponding contact angle 9 P does not 
vary along that contact line. 



(i?o) and of the particle (a) are considered to be suffi- 
ciently large so that in Eq. ([3]) we can neglect line tension 
contributions and Eqs. |T]) and @ do indeed hold [4C| . 

The condition of mechanical equilibrium of the droplet 
in the presence of an external force requires a mecha- 
nism fixing the overall lateral position of the droplet on 
the substrate. The substrate provides a counterbalanc- 
ing force in the vertical direction, but not in the lat- 
eral direction, and therefore the droplet would start to 
move under the action of the external force. There is 
motionless equilibrium only if cither a body-force fixes 
the center of mass of the droplet or, alternatively, the 
contact line at the substrate is pinned. Both situations 



can occur in experiments. For example, in the presence 
of gravity a small tilt of the substrate could provide a lat- 
eral body-force counterbalancing the lateral component 
of fe r . (However, in this case also the weight of the fluid 
must be taken into account.) On the other hand, pin- 
ning of the contact line can occur due to heterogeneities 
of the substrate. To some extent such a pinning is always 
present in actual systems because the substrate is never 
perfectly smooth. Therefore, in practical terms, we ex- 
pect this latter case to be the more relevant one. These 
two different conditions for mechanical equilibrium im- 
pose distinct additional constraints onto the equilibrium 
shape of the droplet and thus onto the minimization of 
the free energy functional in Eq. (0} . We shall use an in- 
dex a in order to distinguish between the cases of a free 
(a = A) or a pinned (er = B) contact line, respectively In 
the limit of small deformations the forthcoming detailed 
analysis will reveal that those two cases correspond to 
Robin and Dirichlct boundary conditions, respectively. 

A fixed lateral position xcm of the center of mass 
(CM) can be achieved by adding a term —Jcm(xcm — 
xcM,ref) to the free energy functional F, where Xcu,ref 
is determined by the reference configuration and —fcM 
is a Lagrange multiplier. Choosing the minus sign in the 
prefactor —fcM will enable us to identify /cm with the 
force acting at the center of mass and fixing it. We in- 
troduce the free energy F\ as a constrained minimum of 
F, with a kept constant and with a fixed center of mass: 

Pl( a 'i fi@o,9 P ,a, Ro, A, fcAi) = 

= min [F - f C Al(xcM ~ XCM.ref)] , (4) 
{r(n)},h 

where the contact angles do and 9 P are given by Eqs. (JTJ) 
and ©, respectively, and Rq = Ro(9 07 9 p ,a 7 Vi) is the 
radius of the droplet in the reference configuration. In 
the following, if not indicated otherwise, we shall sup- 
press the explicit dependence on the set {/, 9q, 9 p , a, Rq} 
of independent system parameters. The Lagrange multi- 
pliers A and fcM can be determined as functions of these 
independent system parameters from the conditions 



V(a; A, Jcm) 
xcM(a;\fcM) 



XCM,ref ■ 



(5) 
(6) 



This renders A = A(a; f,6o,d p ,a, Rq) and 
fcu(oc; f,0o,Op,a,Ro) which upon inser- 

tion into FX(a;f,6 o ,0p,a,Ro,\,fcM) yields 
F A (a;f,6 o ,0 p ,a,Ro). 

On the other hand, the condition of a pinned contact 
line imposes only a geometric constraint onto the mini- 
mization of the free energy, which we symbolically indi- 
cate as 



FZ(a) = min' J", 
{r(Q)},h 



(7) 



a circle corresponding to the one of the reference config- 
uration (Fig. HJa)). Inserting A and Jcm from Eqs. ([5]) 
and ([HI) into F B yields Fb(oi; /, 9 , 9 p , a, Rq). 

In both cases one can split the free energy into two 
parts: 



F a (a) = F aQ + AF (7 (a), 



(8) 



where o = A,B. F a0 :— F a (a = 0) is the free energy 
of the droplet with the adsorbed particle in the axially 
symmetric position a = and AF a := F a — F g q is the 
excess free energy depending on the angular deviation of 
the particle position from the symmetry axis. The latter 
quantity will be the main focus of the following analysis. 



where min' indicates minimizing only over those configu- 
rations for which the contact line at the substrate lies at 



B. Parameterization in terms of spherical 
coordinates and effective description of the particle 

If without a particle the original fluid interface is 
flat, the capillary deformation of the interface around 
an added particle can be described in terms of capillary 
multipolcs, uniquely determined by the deformation of 
the fluid interface outside the particle. Effectively, the 
interface with a particle can be replaced by the interface 
without the particle, but with an infinite set of point- 
multipoles, all centered at a single point inside the region 
originally occupied by the particle [41(. In the following 
analysis we apply a similar approach for the case of a 
spherical interface. The subsequent comparison with full 
numerical results will show that in the model studied 
here the effect of the particle onto the interface can be 
very well approximated by a capillary monopole. Never- 
theless, for the present curved interface, we do not claim 
that there exists an equivalent set of point-multipoles po- 
sitioned at a single point, analogous to the case of a flat 
interface, which is determined uniquely by the deforma- 
tion of the interface outside the particle. Actually, in 
general, owing to the volume constraint, this might not 
be the case [42( . Instead, we assume that the particle can 
be replaced by a set of capillary charges distributed over 
the additional virtual piece of liquid-gas interface inside 
the particle. We introduce a pressure field II(fi) defined 
at the virtual piece of the interface described by the solid 
angle Af2. (In the reference configuration, this piece of 
interface, defined at AQ re f, together with the remaining 
part of the interface, forms a perfect cap of a sphere.) 
Within this approach the contact line at the particle sur- 
face is virtual and enters only through II(f2); therefore, 
in the following, we shall use the notion "contact line" 
exclusively for the contact line at the substrate. In po- 
lar coordinates, one can express the free energy func- 
tional in Eq. @ as a functional of the radial deformation 



u(n) = r(0) - R : 

T[{u(Q)}] = 7 f dSl [s(u, V a u) - R 2 Q ] 

+ 7i?o ( / dCl - I dfl) 

\Ja c \n Jn Q \n c J 

- 2£^ | "^ [ (i?Q + UcW) a sin 2 cW _ ^2 sin 2 ^ 

-i / don(r!)[(i?o + «) 3 -^o] 
+ |[/ dn-/ do] [(i? s (0)) 3 -i?g] 

■J yin c \o Jn \n c J 

- XSV, (9) 

with dfi = d^<i0 sin and 



s(u, V a u) = (R + u) 2 ^l + (V a u)y(R Q + u)2, (10) 
where 



V a :- e e d e 



sini 



A 



(11) 



is the dimensionless angular gradient on the unit sphere. 
We distinguish two integration domains: 

n Q = {(0, 0) G R 2 | G [0, 2i)Aie [0, O ]}, (12) 

which corresponds to the reference droplet shape, and 51 c 
in which the actual droplet interface u(Cl) is dchned: 

n c = {(0, (f>) G M 2 I <f> G [0, 2tt) A G [0, C ($]}. (13) 

fi c differs from fig only by the maximal polar angle 
0c = 6c{4>)i which determines the shape of the contact 
line. Equivalently, the shape of the contact line is also de- 
scribed by the deformation u c = u c (<f>) = u(9 = c (<j>), 4>)- 
In the case of a pinned contact line one has £7 C = f2 and 
u c = 0, but in the case of a free contact line in gen- 
eral Vi c 7^ fly. We note that the parameterization in 
terms of spherical coordinates remains valid for 8 a > n/2 
(in this case the center of the droplet lies above the 
substrate). The first two terms on the rhs of Eq. © 
emerge from 7[/q dO s — L dQ, s re f] with s re f = i? 2 , and 
thus they represent the changes in the surface energy 
of the liquid-gas interface. The third term represents 
the changes in the liquid-substrate surface energy. The 
fourth term represents the work done by the external 
pressure II(J7) in displacing the interface and the last 
three terms correspond to the volume conservation. The 
last but one term corrects the previous one with a liq- 
uid volume wedged between the substrate surface and 
the surface of the cap of the reference sphere inside the 
domains Q c \ fig and f^o \ Q c (in the former case the 
reference surface is extended into the domain Q c \ Oq)- 



Accordingly, R s (9) = i? cos#o/cos# expresses the dis- 
tance between points on the substrate surface and the 
origin in terms of spherical coordinates. The last term 
represents a correction 5V corresponding to the change 
of the volume (due to an interface displacement u) of 
the virtual liquid domain added inside the particle (see 
Fig. [5]). Comparing the expressions in Eqs. © and ^ 
yields an implicit definition of the effective pressure field 
II, in the sense that II(J7) must be chosen such that the 
following equation is fulfilled: 



A<> 



cmiI(ty[(Ro+u) 3 -Rl] = fh+-/SS pl cos6 p 



(14) 




FIG. 2: Schematic cross-sectional representation of the virtual 
piece of the liquid-gas interface (dashed lines) added inside 
the particle on which Il(fi) is defined (e„ is a unit vector 
normal to that interface). The quantity SSi g (see main text) 
is the change in the surface area of that piece relative to the 
reference configuration (u — 0), whereas SV is the volume 
change of the corresponding virtual liquid domain inside the 
particle (hatched regions). The displacement of the center 
of the particle is denoted as h. The radial displacement of 
the interface relative to the reference configuration (Ro) is 
denoted as u. 



where 5S v i = S p i — S p ^ re f; SSi g is the change of the area 
of the virtual liquid-gas interface relative to the refer- 
ence configuration (see Fig. [2]). Equation (JT4J) does not 
uniquely determine 11(17) but in the limit of small parti- 
cles, as will be studied in Subscc. III.D, it provides a con- 
dition sufficient for calculating the shape of the droplet. 



The shift of the center of mass is given by 



Therefore, it is reasonable to choose 



XCM — XCM,ref 



1 



Ra+u 



— I / dft I dr — I d&l I dr r 3 sin 6> cos < 

'l \JQ C Jr b {0) JQo JRs(B) J 
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Wi 
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— / dn[(R + uf - Rl] sin 9 cos 

dfl) [(R s (6)) 4 -R*]sm9cos4 
- 6x, (15) 



<m 



n \n c 



where the second term in the final expression represents 
the contribution to the shift due to the deformation of 
the contact line and Sx is the contribution to the shift of 
the center of mass of the whole droplet due to the shift 
of the virtual liquid domain assigned to the inside of the 
particle. R s (9) has the same meaning as in Eq. (|9|)- 



III. PERTURBATION THEORY FOR LARGE 
DROPLETS 

A. Systematic expansion, variation, boundary 
conditions, and force balance 

It is the aim of this section to base our analytical ap- 
proach on a systematic expansion in terms of a small pa- 
rameter e, such that the deformations of the droplet are 
small, too. We first note that strongly deformed droplets 
are usually metastable or unstable (see also Subsec.V.A), 
so that a necessary condition for the deformations being 
small is the stability of the droplet shape. In the case of a 
spherical droplet this means that the shift in the internal 
pressure due to the force / acting on the interface must 
be much smaller than the Laplace pressure 2 / y/Ro- Up 
to a geometrical factor the shift is of the order of f /Rq, 
so that in order to avoid such instabilities one must have 



l/l « 7^0, 



(16) 



which means that the external force must be much 
smaller than the maximal capillary force acting on the 
contact line at the substrate which is of the order of jRq. 
On the other hand the same reasoning applies to the con- 
tact line at the particle implying that the force / should 
not be stronger than the maximal capillary force acting 
on the particle which is of the order of 7a. One can argue 
that even if the deformation of the interface is large (in 
a not yet precisely defined sense) near the particle, per- 
turbation theory could still be valid far away from the 
particle. Thus, for the perturbation theory to be valid, 
instead of / -C ja one only needs 



l/l<7«- 



(17) 



ML 



(18) 



as a small dimensionlcss parameter, assuming that the 
condition in Eq. (|T7|) is fulfilled, too. We note that for 
large droplets, i.e., for 



a<i?o, 



(19) 



the condition e <C 1 is automatically satisfied. 

In the next step, we expand u, A, II, feu, and, for 
later purposes, h in terms of the small parameter e: 



u(a) = R (ev(a) + o(e)), 

A=-^(2 + e M + o(£)), 



R 



7 



n(n) = ^-(c7r(n) + o(e)), 

/CM = jRo{eQcM + o(e)), 
h = R Q (eh + o{e)), 



(20) 

(21) 

(22) 

(23) 
(24) 



where A has been expanded around the Laplace pressure 
2j/Rq of a spherical droplet; v, /i, w, Qcm, and h are 
dimensionlcss. 

In the following we shall keep the radius i?o of the 
droplet fixed and therefore it is convenient to divide the 
energy by jRq. In zeroth order in e the droplet is un- 
deformed (u = 0) and the free energy functional T^ as 
defined in Eq. ([9]) equals zero. Similarly, the free energy 
J 7 ^ in first order in e also vanishes, which reflects the 
fact that we have chosen the equilibrium configuration as 
the reference one (in equilibrium, by definition, pertur- 
bations do not give rise to linear contributions to the free 
energy). Therefore the leading contribution is of second 
order in e and the free energy functional reads 



T 
7^g 



r/O 



-(Vav) 2 -v 2 -(n(n) + Li)i 



■ COS #o 



2tt 



■ H eo ) 2 + O(e 3 ), (25) 



where tt(D, (fc Ail) = and the boundary term (second 
term) has been obtained by expanding up to second or- 
der in e those terms in the free energy functional in Eq. 
©, which depend on the deformations of the contact 
line. (For 9 c (<f)) < 6q, the deformation v(8) is smoothly 
extended to 9 = 9 , so that in Eq. (|25j) v(8) is defined 
within the whole angular domain Q,q ; this is justified be- 
cause the corresponding difference in the free energy is 



of higher order 0(e 3 ).) The variation of T renders 
[T\p + Sv] - F[v]) = 



jR 2 o 



■ e 2 J <m [-V 2 v -2v- tt(O) - n]6v 
Jn 

+ e 2 / dnW a ■ (Sv\/ a v) 

Jn a 
p2lT 

-e 2 cos8 Q d<j>v\g 5v\ 9o +0({Sv) 2 ,e 3 ), (26) 
Jo 

whereas the variation of the additional term responsible 
for fixing the center of mass (Eq. (fT5|l) yields 



fcM 

iRl 



(XCM — %CM,ref)\ 



v-\-8v 



e 2 3QcM , I dnSvsmdcoscfi + OaSv) 2 ^ 3 ), 



47r/o(6»o) Ja 



(27) 



where we have neglected the second term in Eq. (|15|). 
which is of the order 0(e 2 ) and thus in Eq. ([27| it con- 
tributes only to the order 0(e 3 ), and Sx which con- 
tributes only to the order 0(e 2 ) x O(a/R ) 3 . The func- 
tion /o(#o) expresses the volume of a unit spherical cap 
characterized by the polar angle 8q as a fraction of the 
volume 47r/3 of a unit sphere: 



fo(x) := (2 + coscc)sin (x/2). 



(28) 



The expressions in Eqs. ([2"B)) and (f2"7| lead to the Euler- 
Lagrange equation for a droplet with a fixed center of 
mass in the form 



(v 2 + 2)«(n) = 7r(n) + 7r CM (f2) + /i, (29) 



where 



ttcm(^ = sm6>cos( 

47r/o(c'o) 



(30) 



is an effective pressure fixing the center of mass of the 
droplet. In the case that the center of mass is not fixed, 
the corresponding Euler-Lagrange equation has the same 
form as in Eq. ([2"9")l but with tt C m = 0. Applying Gaufi' 
theorem (on a unit sphere) to the second term in Eq. (|26|) 
gives the &oundary contribution 



1 



jR 2(?lv + Sv}-?[v}) bc 



2n 



d4> (sin 8 d e v\g - cos8vv\g )&v\g 01 ( 31 ) 



vanishing if 5v\$ = 0, which corresponds to a pinned 
contact line (model B). If the contact line is free (model 
A) one has 5v\e ^ and, instead, the expression in 
brackets of the integrand in Eq. pip must vanish. These 



two different conditions correspond to Dirichlet (B) and 
Robin (A) boundary conditions, respectively: 



A (free) : 
B (pinned) 



sm9 d e v\g - cos v\ 8o = 0, (32) 

v\g = 0. (33) 



Equation (|32|) is equivalent to the condition that the 
contact angle is equal to the Young angle, which can 
be seen by the following reasoning. In the limit of 
small deformations the normal e„ to the droplet sur- 
face in terms of spherical coordinates can be expressed as 
e„ = e r — eV a v + 0(e 2 ), so that the contact angle 9(4>) 
at the substrate is given by 

cos6>(^6) = e z -e n \g a = cos8 c + esm6 c d g v\g c +0(e 2 ) = 
= cos8 + e[sm9 d e v\g -v\g cos8 Q ] + 0(e 2 ), (34) 

where in the third equality we used cos8 c = cos#o(l — 
ev\g ) + 0(e 2 ), which follows from the analysis of small 
perturbations of a spherical cap. Thus, from Eqs. (|32|) 
and (f3"4"| we recover Young's law in the form 6{<f) = 9q. 
We note that whereas the condition in Eq. (|33|) corre- 
sponds to the shape of the droplet at the actual contact 
line, because for a pinned contact line C (4>) = 8q, the 
condition in Eq. (f3"2"j) does not apply directly at the con- 
tact line dfl c but at the boundary <9f2o of the reference 
integration domain. The actual shape of the contact line 
in this case can be obtained by a linear extrapolation of 
the interface profile from 6® to 8 C . 

The volume constraint in Eq. ([5]) in first order in e 
reads 



dQv(n) =0, 



(35) 



n„ 



whereas in Eq. ^ the constraint of fixed center of mass, 
due to Eqs. (fl~5|) and (|30p . can be written, up to second 
order in e, as 



dTlircM(to)v(n) = 0. 



(36) 



In Eqs. (|3"5"|) and (|36[) we have neglected contributions 
from SV (Eq. (J9])) and 8x (Eq. (fTSjl 'l. which are both of 
the order O(a/R ) 3 . 

In the following our goal is to express the pressure shift 
/i and the balance of forces acting on the droplet in terms 
of 7r and 9q. Integrating both sides of Eq. (f2"9"|) over Oo, 
with J n dQ = 2tt(1 — cos8 ), and applying Gaufi' theo- 
rem and subsequently the condition of constant volume 
(Eq. dMD) yields 



27r(l-cos6» ) 



dfl tt (il.) + sin 8 o / difidgv 



n„ 



I On 



(37) 

where we have used the fact that L, dft ttqai = 0. Equa- 
tion (|37[) can be interpreted as a hydrostatic balance in 
that the internal pressure (lhs) is equal to the external 



pressure exerted by the forces acting at the droplet (rhs) , 
which in particular depend on the boundary conditions. 
The governing equation (|29|) can also be used to derive 
the force balance on the droplet. First, multiplying both 
sides by the radial vector e r (Q) and integrating over £Iq 
we obtain 

- / dne r (n)(V 2 a + 2)v = 
Jn 

= / d£lTr(Q)e r (Q) + QcM e x + ^7rsin 2 9 e z . (38) 



The first term on the rhs of Eq. (|38| represents, up to 
first order in e, the total external force acting on the 
droplet. The surface integral on the lhs, after integrating 
by parts, using the fact that (V 2 + 2)e r (Q) = 0, which 
follows from the identity (V 2 + 2)Yi m = 0, and finally 
applying Gaufl' theorem can be transformed into a line 
integral, 



/ dne r (n)(\7 2 a + 2)v 



/■27T 

= -sin9 d(j)(e r d e v -vd s e r )\ 0Q . (39) 
Jo 



Using dge r = eg and taking the limit a/ Rq — > Eq. 
can be expressed in terms of Cartesian coordinates: 



x : 



sin 6*o 



cos4>(sm9 dgv\g -v\g cos9 ) 

= / d£L tt(Q) sin 9 cos cf> + Qcm, 
Ju 



(40) 



y : — sin^o / d(f> s'm<f)(sm9odgv\g — v\g cos9o) 



z:—sin9 / dcj) ( cos9 dgv\g + v\g s'm9 ) 



/ d^l tt(^1) cos 9 + fj,n sin 2 
Jn 



= 0, 

(41) 

(42) 



70- 



In the case of a free contact line (see Eq. (f3"2j)') the left 
hand sides of Eqs. (|40| and (|4Tj) vanish and the force 
balance in the x-dircction reduces to 



Qcm = / dfln(fl) sin#coS( 
Jn 



(43) 



If 7r(0) corresponds to a pointlike particle (see, c.f., Eq. 
dnOJ)), so that tt(O) = ^(O,0i) = qS(9 - 9 l )8{(j) - 
(f>i)/sin9, where f2i = (9\ = a,<j>i = 0) is the direc- 
tion along which the external force pulls (q = +1) or 
pushes (q = —1) the particle, due to fcM = ejRoQcM = 
|/| Qcm (sec Eqs. (HHJ) and (E3D) Eq. I® leads to f CM = 
—/sin a. Accordingly, the Lagrange multiplier /cm (Eq. 



(|4|)) can indeed be interpreted as a force counterbalancing 
the x-component of the external force, equal to / sin a, 
and thus fixing the center of mass of the droplet. In the 
same limiting case the first term on the rhs of Eq. ([55)1 
reduces to J n d£lir(£l)e r (£l) = qe r (fli) corresponding to 
the external force on the droplet. In the case of a pinned 
contact line and with a free center of mass (i.e., fcM = 
so that Qcm = 0) the force balance in the ^-direction 
(Eq. (HOD) reads 

— sin 2 6*o / dcj) cos (j)dgv\g Q = / dSl n (SI) sin 9 cos 4>. 
Jo Jn 

(44) 
Thus, in this case, the x-component of the external force 
(rhs) is counterbalanced by the capillary force due to the 
deformation of the droplet at the contact line (lhs) . 

For model A the lhs of Eq. (|4i"j) vanishes due to Eq. 
(|3"!?|) . but it also vanishes for model B because dgv\g 
and v\g are symmetric and sin</> is antisymmetric with 
respect to the xz-pl&nc. 

The difference between these two models is that for 
model A the local contact angle is the Young angle (see 
Eqs. (|34[) and (f3"2| ). which guarantees that each piece of 
the contact line is in mechanical equilibrium, whereas in 
the case of model B this is not the case and the total 
capillary force on the contact line counterbalances the 
x-component of the external force (Eq. (|4"4"|0 . 

The boundary condition in Eq. (J35J) allows one to ex- 
press dgv\g in terms of v\g . Inserting this into Eqs. 
(|3"T|) and (|4*2")l renders two equations for [i and J^ d(f> v\g , 
which can be solved to yield 



A : /i 



B : n 



1 



47r/ (6»o) Jn 

1 



<iil 7r(il)(l — cos 6*o cos( 



7r(l — cos#o) 



driTT(n)(cos9 — cos9q). 



(45) 

)• 
(46) 



Equation (|4^|) follows from inserting the boundary con- 
dition in Eq. (|33|) into Eq. (J42|) : this renders L dcf>dgv\g 
in terms of /i, which can be inserted into Eq. (|37[) 
yielding an equation for /i. In the case of a point- 
like particle (see, c.f., Eq. (J60j) ) , Eqs. gS]) and (gSl) re- 
duce to fi = —q(l — cos9 cosa)/(4irf (9 )) and fj, = 
—q(cosa — cosflo)/( 7r (l — cos^o) 2 ), respectively, so that 
due to a < 9q (see Fig. HJa)) the internal uniform pres- 
sure shift fi has always the sign — q, i.e., the opposite one 
to /. This means that the pressure change counteracts 
the action of the external force, which can be interpreted 
as the realization of Le Chatelier's principle for this par- 
ticular system. 



B. Green's functions 

In this subsection we present a formal solution of Eq. 
(|2"9")l by means of Green's functions G,j(£l,n') which are 



the radial deformations t>(fi) due to a pointlikc pertur- 
bation 7r(fi) = <5(fi,fi') = 6(0 - 6')5{<f> - <f>')/ sin 6 act- 
ing in the direction fi' — ($', (/)'). In this respect we 
point out that in the case of arbitrary fi' (i.e., a direc- 
tion pointing out of the a;z-plane) the term —Jcm{xcm~ 
xcM,ref) in the free energy functional in Eq. (j4]) has to 
be replaced by -fcM, x (xCM - xcM.ref) ~ fcM,y(yCM ~ 
VCM,ref) where fcM,x and fcM,y are the correspond- 
ing Lagrange multipliers. Accordingly, in analogy to Eq. 
(f3"0")) . in this case one has ttcm(^) = 3 sin 0[Qcm,x cos(j)+ 
Qcm, v sin 4>}/(4Trfo(9o)) with Q C m,x ■= fcM,x/\f\ = 
Qcm given by Eq. (gSJ) and Qcm, v ■= fcM, y /\f\ = 
— Jq dfi 7r(0) sin sin 0. This means that Green's func- 
tions fulfill Eq. (JS5]) with 7r CA/ (£l) = 3 sin 6 sin 6' cos(0 - 
(/>')/ (4irfo(0o)) for model A, ttcm = for model B, and 
/i given by Eqs. (|4"5"j) and (|46|) . Accordingly Green's func- 
tions fulfill 

-(yl + 2)G A (n, n') = d(n,n') 

3 ./,/•/, f i ,,s l-cos6»ocos6»' 

sin sin 6* cos(</> — cp ) — 



^fo(0 ) 



47r/ o (0o) 



(47) 



for model A and 



-^ + 2)05(0,0'; 



,, COS0' — COS0n . . 

= ( Jfi,fi' -_ — § 48 

7T(1 — COS0 O ) 



for model B, where in both cases fi, fi' € fio. We note 
that for arbitrary 9q the Green's functions are not sym- 
metric, i.e., G CT (fi, fi') ^ G CT (fi' , fi). This means that the 
deformation at fi due to a perturbation at fi' in general 
differs from the deformation at fi' due to a perturbation 
at fi. This can be understood by the fact that the shift p 
of the internal pressure depends on where the perturba- 
tion is applied. The only exception is the case 6$ = 7r/2 
with a free contact line, for which the Green's function 
Ga is fully symmetric. In this case cos 9q = so that the 
last term on the rhs of Eq. (|4T|) is a constant, whereas 
the last but one term is explicitly symmetric. 

The deformation of the droplet for an arbitrary pres- 
sure field 7r(fi) can be expressed in terms of the Green's 
functions: 



v a (Cl)= dil' G a (il,il')n(n'). (49) 



By using Eqs. (|47)) and ((48J) one can check that this ex- 
pression for v a fulfills the Young-Laplace equation (|2U1) . 
For 7r(fi) = (5(fi, fi') the volume constraint in Eq. ([33)) 
yields 

J dfiG CT (fi,fi') = 0, (50) 

Jn 

whereas the boundary conditions in Eqs. (|3"2"j) and (f3"3")l 



can be expressed as 

sin0 o 9eGA(fi,fi')|nean o - cos0 o GA(fi,fi')|nean o =0, 

(51) 

G B (fi,fi')| ne an =0. 

(52) 



C. Free energy 

It can be shown that for the boundary conditions of 
either a pinned or a free contact line the total boundary 
contribution to the free energy vanishes (see Appendix 
A). In analogy to electrostatics (see, for example, p. 46 
in Ref. [43|) we obtain the following expression (derived 
in Appendix A) for the free energy of a sessile droplet 
subjected to an arbitrary external pressure 7r(fi): 



F„ 



dfi7r(fi)w(fi) 



27 JO, 
- — [ dfl I dfi / 7r(fi)G (T (fi,fi')7r(fi'). (53) 



where for the second equality we have used the general 
form of the solution given in Eq. (|49|) . 



D. Limit of small particles 

Taking the limit a — > allows one to provide a relation 
between the unknown pressure field 7r and the system pa- 
rameters introduced in Sec. II. In this context, Eq. (|14|) 
provides a constraint on tt expressed in terms of an un- 
known deformation field, encoded by the quantities h, 
SSpi, and 6Si g . However, this equation simplifies con- 
siderably in the limit of small particles. In order to see 
this, we consider Eqs. (|57|) and (|4"8"j) in the limit fi — > fi' 
in which the spherical reference interface becomes locally 
flat in the neighborhood of fi and fi'. First, we rotate the 
coordinate frame such that the z-axis points in an arbi- 
trarily chosen direction in the neighborhood of fi and fi'. 
(In the following we refer to this new axis as the z-axis.) 
Next, we project the points fi and fi' on the unit sphere 
along the z-axis onto the plane tangent to the reference 
spherical interface at the point of intersection of this in- 
terface with the z-axis. The mapping (0,0) i-> (p,4>), 
where (0, <fi) denote the spherical coordinates associated 
with the z-axis and (p, 4>) are polar coordinates in the 
tangent plane, is given by 



p = B{) sin 1 



(54) 



In order to keep the notation simple, in the following we 
skip the tilde so that (p, (j>) = (p, 4>). Under the transfor- 
mation of coordinates the Laplace operator transforms 
according to the chain rule whereas the delta function 
transforms as: S((f> — <f>')6(8 — 0')/sin0 = RqS(4> — 
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0') 6(0(p) - e 1 ((/))/ p = R \d P /de\ s^ - </>') s( P - p')/ P . 

Accordingly, Eqs. (|47p and (|48p can be rewritten in terms 
of a single equation (for details see Ref. |44j): 

- {R 2 Vf + 2 - 9 P p 2 d p )G a {x, x 1 ) = 



R ( 



R l 



p 2 6(x-x')+A a (x,x'), (55) 



where x denotes the Cartesian coordinates within the 
tangent plane, Vm is the Laplace operator on the tan- 
gent plane, and G a (x, x 1 ) :~ G a (Q(x), Q'(x')); the func- 
tion A a (x,x') := A a (n(x),H,'(x')) stands for all non- 
singular terms on the right hand sides of Eqs. (|47|) and 
(|4"8)l . where x(fl) = (x = Rq sin 6* cos 0, y — Rq sin 9 sin (f>) 
describes the rotation O i — >- O, which might be expressed 
in terms of the corresponding Euler angles, with the sub- 
sequent projection onto the tangent plane fi i— > x = x 
(see also Eq. (|5"4")l ). In the limit R — >• oo Eq. (|55|) reduces 
to the usual two-dimensional Green's equation 



V\G a {x,x') =6(x-x'), 



(56) 



with the solution G a (x,x') = G(x,x') = G(\x — x'\) = 
— (1 / (2ir)) \n(\x — x'\) which depends neither on a nor 
separately on x' , or equivalcntly Q', i.e., where on the 
droplet the perturbation is applied. 

Therefore the deformation v a diverges logarithmically 
in the neighborhood of a pointlike perturbation at ft' , 
independently of the boundary conditions far away from 
Q,'. By taking 7r(0) = d(Sl,Q') in Eq. (gSJ one obtains 



O — vO' STT 



n->n' 



2tt 



(57) 



where 9 is the angle between the unit vectors pointing 
into the directions f2 and £1'. Thus, with the radius a of 
the particle acting as a natural cutoff, in leading order 
in a/Ro the dimensionless displacement h of the particle 
(see Eq. (|2"4"j)) reads 



2tt 



0(1), 



(58) 



where q = //|/| = sgn(/). Accordingly Eq. (fT4|) in lead- 
ing order in e can be written as (see Eqs. (f20|) . (|22|) . (|24| . 
and (fT8]>) 



/ rffi 7r(fi)v(fi) = — In | — ) + 0(1 
•/n 2tt V a / 



(59) 



The correction term in Eq. (|59|). which stems from the 
last two terms on the rhs of Eq. (Q3)), is proportional to 
7 2 a 2 // 2 which, in agreement with Eq. (JTTJ), is of order 1 
relative to the leading term which diverges ~ hx(Ro/a). 
Moreover, assuming that the virtual piece of the inter- 
face defined at AJ7 is bounded to lie inside the particle, 
one has u(0)|ao — h + O(a/Ro) (see Fig. [2]) so that 
v(CI)\aq = h + O(a/Ro). According to Eqs. ((551) and 
(1551 this leads to 



dflir{n) = q + O{l/ln(R /a)). 



(60) 



Equations ((59)) and (|60|) are satisfied by tt(O) = 
g(5(fi,fii) with O x = (0 = a,<£ = 0) (see Fig. p. Thus, 
in leading order in a/Ro the integrated amplitude of the 
effective pressure is independent of the particle radius. 
In this limit the droplet radius Rq is the only remaining 
length scale. 

According to the above considerations, for — > Q,' 
Green's functions G a (n, f2') in Eq. ([SU)) contain a loga- 
rithmically diverging part G(il,fl') ~ ln0. Therefore it 
is useful to consider this singular contribution separately 
by writing G a = G + G a . reg , such that the regular part 
G a ,reg = G a — G does not diverge for ft — > ft'. (As de- 
scribed in the subsequent section in special cases G can 
be identified as Green's function for a free droplet (see, 
c.f., Eq. @5)).) 

With 7r(0) = q S(il, Oi) this leads to the following ex- 
pressions for the free energy (Eq. (fB"5|) ): 



F, 



o-O 



and 



AF„ 



f 2 
~ [G(0 = a/Ro) + G^reg(0, 0) + 0(1)] (61) 



r- 



■ [G a , reg <fl u fii) - G CT , res (0, 0) + O(o/iio)] 

(62) 
where G(51 = fii, O' = S7i) for Oi = (the apex position) 
has been regularized as G(9 = a/Ro). The leading term 
of the free energy F a o stems from the singular part of the 
kernel in (Eq. ([55)) ) and thus, in analogy to electrostat- 
ics, it can be interpreted as a self-energy of the particle. 
However, unlike the self-energy of a pointlike charge in 
electrostatics, this quantity does not diverge because the 
singularity is suppressed by the prefactor f 2 Eq. (|6l)) the 
upper bound of which, as already noted, is proportional 
to 7 2 a 2 (Eq. (JTTJ)). Thus for a -> with //(7a) kept con- 
stant the self-energy F a o vanishes as a 2 ln(Ro/a). In the 
following section we shall rather focus on the remaining 
part of the free energy AF a which determines the angular 
free energy landscape of the particle. 

Comparing the first equality in Eq. ((53)) with Eq. ((59)) 
one can see that a correction of order O(l) due to the 
surface energy of the particle (see Eq. (fT4")l ) also con- 
tributes to the free energy F a . As indicated in Eqs. (fBT)) 
and (f62)) . for "fF a o/ f 2 this correction contributes to the 
order O(l) but for jAF a /f 2 it contributes only to the 
order O(a/Ro) or higher. This is the case because the 
changes of the surface energy of the particle with respect 
to a occur only due to the presence of the substrate (the 
free energy involving a free droplet would not depend on 
a) which is located at a distance O(Ro) from the particle 
and thus gives rise to the above factor O(a/R ). 



IV. METHOD OF IMAGES FOR o = tt/2 

For the special case 9 = n/2 one can construct the 
Green's functions G a by using the method of images. 
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First, in the reference configuration, the substrate is re- 
placed by a virtual mirror image (with respect to the 
surface plane of the substrate) of the actual reference 
droplet. We obtain a full and free sphere composed 
of the actual upper hemisphere and the virtual lower 
hemisphere. Next, the deformation of the actual droplet 
subjected to a pointlikc force can be constructed as the 
upper part of the full droplet subjected to this point- 
like force plus additional virtual pointlike forces, subse- 
quently called images, placed at the virtual lower hemi- 
sphere. The distribution of these images has to be chosen 
such that the following conditions are satisfied: 

(i) the shape of the droplet, given at a direction ft by 
G CT (ft,ft'), obeys the boundary conditions in Eqs. 
([STj) or ([Hll), depending on the model; 

(ii) the volume of the actual sessile droplet is conserved 

(Eq. (EOD); 

(iii) the total force on the full droplet vanishes in the 
case of a pinned contact line; in the case of a free 
contact line the force balance is automatically sat- 
isfied by fixing the center of mass. 

Condition (iii) represents a sufficient condition for me- 
chanical equilibrium of the sessile droplet, because it im- 
plies that every piece of the full spherical interface is in 
equilibrium. In fact, conditions (ii) and (iii) are auto- 
matically satisfied for the Green's functions which obey 
Eqs. (|47|) and (|48|) . because those equations have been ac- 
tually obtained under the conditions of force balance and 
volume constraint. However, the physical assumptions as 
formulated in (ii) and (iii) can provide a guide for find- 
ing the distribution of images. The boundary condition 
expressed in (i) has to be imposed separately. 

According to the above reasoning, the Green's func- 
tions for #o = 7r/2 can be written in the following general 
form: 

G CT (ft, ft') = G(n, n')+J2 Q*iG(n, n ai )+G„, corr (n, ft'), 

(63) 
where G(ft,ft') = G(ft',ft) (see, c.f., Eq. flB5J|) is the 
Green's function for a free droplet (see, c.f., Eq. (|6"4"]l). 
Both the amplitudes Q„i and the directions ft,^ of the 
images, as well as a possible correction G a ^ corr must be 
chosen such that the conditions (i) — (iii) and the Green's 
equations (|47|) and (|4"5)) arc satisfied. The problem of a 
single pointlikc force acting at the surface of a free droplet 
with a fixed center of mass has been studied by Morse and 
Wittcn in Rcf. [431, who derived the following equation 
for the free G: 

-(V* + 2)G(ft,ft') = J2 Y l * m (Q)Y lm (Q') = S(Q,Q'), 

l>2,m 

(64) 
where the rhs is a modified delta distribution <j(ft, ft') 
with the components I = and I = 1 projected out, so 
that G(ft,ft') = £ J>a _(l(J + 1) - 2)- 1 17 m (ft)Y hn (ft'). 



This reflects the conditions of incomprcssibility of the 
liquid (f 3 dftG(ft,ft') = where S2 is the full unit 
sphere) and of balance of forces acting on the droplet 
(J S9 dft e r G(ft, ft') = 0). The solution of Eq. flM]) can be 
expressed in closed form as [45j 



G(ft, ft') = G(6) = 



1 

47T 



cos 9 + cos 9 In 



(65) 



where 9 is the angle between the unit vectors pointing 
into the directions ft and ft'. Morse and Witten pointed 
out that the deformation v of a full droplet can be de- 
scribed in terms of G also in the case that the center of 
mass is not fixed but, instead, there are other pointlikc 
forces distributed such that the total force on the droplet 
vanishes. We shall use this latter property of the free 
G for determining the Green's function Gb of a sessile 
droplet with a pinned contact line. 

As we shall show in the following subsections, in the 
case 6*o = tt/2 the method of images can be applied suc- 
cessfully because the reflection at the surface plane of the 
substrate provides a smooth interface between the actual 
and virtual droplets. In the case 9 ^ 7r/2 the inter- 
face has a non-physical cusp and the method cannot be 
applied directly. However, we cannot rule out the possi- 
bility that there is a modification of the method which 
can be used successfully also for 9q ^ tt/2 (for the case 
a = sec Appendix C). 



A. Free contact line 

In the case of a free contact line a sessile droplet would 
move along the substrate if exposed to a non-vanishing 
lateral component of the external force applied to the 
particle. Therefore, in order to achieve a motion-free 
equilibrium, for example the lateral position of the center 
of mass has to be fixed. Physically, this compensation 
can be achieved by applying a body force to the droplet. 
Fixing the center of mass is to a certain extent artificial, 
but provides the simplest model for a free contact line. 
(For example, if a particle would hit the droplet and the 
relaxation of the position of the center of mass would 
last much longer then the relaxation of capillary waves 
generated by the impact, then the inertia of the droplet 
would effectively fix the center of mass. Alternatively, 
in the presence of gravity and depending on a a small 
tilt of the substrate could provide a lateral body force 
counterbalancing the lateral component of fe r .) 

For 9 = 7r/2 the boundary condition in Eq. (|5"Tj) takes 
the simple form 



3flGA(fi,n')|neen„=0, 



(66) 



which can be identified with the Neumann boundary 
condition. Equation (|66j) and the volume constraint 

L^dff) L d6* sin 0Ga (ft, ft') = can be satisfied by the 
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Green's function given in Eq. (|63l) with Qa\ = 1 and 
S1a\ = Zf2' where Z means the reflection with respect to 
the plane z = (see Fig. [3|): 



G A (n, n') = G(fi, fi') + G(n, zn') 



(67) 



If the /u// drop (union of the actual and of the virtual 
drop) is exposed to two point forces 7r(f2) = S(Cl,il') + 
S(fl, ZQ') its resulting shape, given by Eqs. (|49)l and (|67| . 
renders the shape of the actual droplet for z > 0. For this 
ansatz, with G(Q, 0') given by Eq. ([65]) . it can be checked 
that OsGa vanishes at the contact line. The correction 
term G A,corr vanishes, because the contributions to the 
volume stemming from the two terms in Eq. (|67j) cancel 
each other, so that the volume constraint is fulfilled for 

*~rA,corr = U. 

In order to discuss the force balance we consider, with- 
out loosing generality, <fi' = 0. In such a case, due to the 
mirror symmetry with respect to the xz-plane and the 
mirror symmetry with respect to the icy-plane, the y- 
component and the z-component of the external force 
on the full drop due to the two point forces vanish, 
whereas the x-component has the strength 2 sin 9' (see 
Fig. [3] with a = 9' and / = 1). Mechanical equilibrium 
in this direction must be restored by fixing the center of 
mass which corresponds to applying the forces of strength 
Qcm = — sin9' (Eq. (|43|) ) both to the actual and to the 
virtual drop. 

Finally, by using Eq. (|64|) one can check explicitly that 
Ga as given by Eq. (|fJT|) satisfies Green's equation for 
9 = tt/2. 



B. Pinned contact line 

In this subsection we consider the case that the con- 
tact line is pinned at the circle corresponding to the ref- 
erence configuration. One can propose that this pinning 
can be accomplished by substrate heterogeneities; but 
this requires a dedicated design in order to maintain the 
circular shape of the contact line. This model has the 
virtue that in contrast to the previous one the balance of 
forces is automatically satisfied without fixing the center 
of mass. 

By choosing Qbi = — 1 positioned at Qbi = Zil' the 
boundary condition in Eq. (|52[) is automatically satisfied. 
However, the corresponding pressure field for pointlike 
forces placed at f2' and at ZQ,' , which corresponds to 
7r(0) = 8(0,, Q') — S(Ct,ZCl'), leads to a non-vanishing 
net force acting on the droplet as a whole, which is di- 
rected vertically and is equal to 2cos#' (see Fig. [4j. In 
order to restore the force balance (without fixing the cen- 
ter of mass) we place a second image Qb2 = 2 cos 0' at 
the bottom ("south pole", 9 = tt) of the virtual droplet, 
so that n B 2 = ^tt (see Fig. [4] with / = 1 and a = 9'). 
However, this second image contributes to a deformation 
at 9 = 9q = tt/2, which violates the boundary condition 
in Eq. (|52[) . This additional contribution, which is equal 




-/sin a 




f 



FIG. 3: In the case of a free contact line the total drop, being 
the union of the actual (upper full line) and virtual (lower 
dashed line) drop, is pulled by / in the upper part and by 
its mirror image in the lower part, which corresponds to the 
Green's function in Eq. (|B7J (with Q.' — Qi = (0' — a, <f>' — 
0)). The lateral force of magnitude —/sin a fixes the center 
of mass ( x ) of each hemisphere separately. 



to 2(cos9')G(fl = (9 = TT/2,<f)),n 7I ) =: -1(9') and thus 
depends neither on Q nor on <fi' , has to be subtracted from 
the Green's function in order to uphold Eq. (f52"j) . Finally, 
the volume constraint is fulfilled by adding a second term 
H(Cl') cos 9, which corresponds to a rigid vertical transla- 
tion (i.e., in z direction) of the droplet and to an effective 
pressure -{V 2 a + 2)H(n')cos9 = -~H(n')(V 2 a + 2)cos6» 
which vanishes independently of H(Q') (which is in ac- 
cordance with a rigid translation). Thus, this term does 
not contribute to the net force on the droplet. Moreover, 
it vanishes for 9 = tt/2, i.e., at the contact line, so that 
it does not violate the boundary condition in Eq. (|52p . 
Accordingly, we are led to the following ansatz for the 
correction term: 

G B , corr (0, Q') = tf(fi') cos (9 + 1(6'), (68) 

so that the ansatz for the total Green's function is 

c? B (fi, n') = G(n, a') - G(n, zn') + 2G(n, n w ) cos e' 

+ H(Q')cos9 + I(9'). (69) 

Inserting the above expression for Gb into the boundary 
condition in Eq. (|52[) renders 



I(x) 



cosx 

Att 



(70) 



The volume constraint in Eq. (|50l) provides an expression 
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-/ / 



-G> 




2/ cos a 



FIG. 4: Distribution of images in the case of a pinned contact 
line, corresponding to the Green's function Gb in Eq. (]69p 
(with Q! — Or = (8' — a, <j> = 0)). The correction Gb,cott 
(Eq. (|68p ) does not give rise to any forces. The center of mass 

is free. 



for F(O') =H{B')\ 



H(x) 



1 

2^ 



cos x In 



1 + cos x 



(71) 



Finally, by using Eq. (|70|). it can be checked explicitly 
that Gs as given by Eq. (|qTt|) satisfies the Green's equa- 
tion (gg]) for O = tt/2; the contribution to Eq. (gSJ due to 
#(fi') (see Eq. ([521)) vanishes because (V^ + 2) cos # = 0. 



C. Free energy 

As discussed in Subsec. III.D the excess free energy 
AF a is well-defined even in the limiting case that 7r(fi) 
is of the form 7r(fi) = <7($(f2,f2i). In this case and for 
6a = it/ 2 one obtains the closed expression (Eqs. (J62J , 
HS3 and dMD) 



AF„(a) 



Z! 

27 



[3(7 («) -,9o 



(72) 



where the functions qa and gs are both regular at a = 
and given cxplicitely by 



9A (a) = G{6 = ir- 2a) 



(73) 



and 



g B (a) = -G{6 = it - 2a) + 2G(d = n- a) cos a 

+ H (a) cos a + I (a), (74) 



with G{6) given by Eq. (|64|) ; the functions /(a) and H(a) 
are given by Eqs. ([70]) and (f7Tj) . respectively. 



V. NUMERICAL RESULTS 

In this section we compare the analytical expressions 
for AFa and AFb, given within perturbation theory by 
Eqs. (J7"2"l) - (f74|) and ([65]) for point force induced deforma- 
tions, with the numerical minimization of AF defined by 
Eqs. I©-© for a deformation induced by a particle of 
finite extent a (Fig. [1]) and for 9$ = tt/2, employing a 
finite element method [46|. The contact line at the par- 
ticle is taken to be free which leads to a fixed contact 
angle 8 p at the particle, in analogy to a free contact line 
at the substrate which fixes the contact angle 8 at the 
substrate. In the case of axially symmetric configura- 
tions we compare the numerical results with the exact 
expressions given in Appendix B. 



Axially symmetric configurations with free 
contact line 



First, we study the axially symmetric configuration 
with the particle positioned at the apex of the droplet. 
For this particular configuration the free energy can be 
calculated analytically solving the full non-linear Young- 
Laplace equation, i.e., without linearizing in terms of e 
(Appendix B). In this case it is more suitable to choose 
the displacement h instead of the force / as the inde- 
pendent variable. The free energy FAo(f) for a = 0, as 
introduced in Eq. ([8]) and in the text preceding Eq. ((5J, 
can be expressed as the Legendre transform of the free 
energy FAo(h) given implicitly by 



FAo(f,0 ,9 p ,a,R ) 



mm[FAo{h,6 , 

h 



9 p ,a,Ro)-fh}. 

(75) 
In the following, when referring to "stability" of the 
branches of the free energy (Fig. [5]) we shall actually 
mean the stability of the configurations with an external 
force / acting on the particle (with / = dFAo/dh, see 
also Appendix C). Another choice for the independent 
variable might be the angular position of the contact line 
at the particle parameterized by the polar angle j3 (see, 
c.f., Fig. \TT\b) in Appendix B). In this case and for a 
free contact line at the substrate the expression for the 
free energy in terms of elliptic functions is derived in Ap- 
pendix B. We have carried out finite clement calculations 
also for Fao(K) so that the comparison with the analytic 
results serves as a welcome test for the performance of 
the numerical code we have to rely on for configurations 
which are not axially symmetric. (Note that for a = 
there is no need to fix the lateral position of the center 
of mass.) In this case the numerical procedure consists 
of employing the finite element method |46[ for the min- 
imization (for details see Subsec. V.B) of the free energy 
functional in Eq. ([3]) with a = 0, / = 0, and for fixed h. 
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We have found very good agreement between the numer- 
ical and analytical results (see Fig. [5]). 

In Fig. [5]we present the results for 6q = 7r/3, 9 p = tt/2, 
and Ro/a f=a 8. One can distinguish seven branches of 
the free energy both for h < and for h > 0. The 
metastable branches denoted as 7 lead to a maximal value 
of Fao on each side of the curve and they correspond 
to configurations with the contact line close to one of 
the poles of the particle. Once the system is prepared 
in such a configuration at the upper part of branch 7 
and the constraint of fixed h is lifted, the particle slides 
down branch 7 by shrinking the contact line while moving 
closer to the reference interface h = 0, ending up in a 
state in which it is fully immersed in the liquid or in the 
gas phase. 

The filled squares in Fig. [5] indicate inflection points 
of the free energy curve. According to the following rea- 
soning, for a large droplet these points are located (for 
6 p = 7r/2) close to f3 = tt/2 — tt/4 = 7r/4, which cor- 
responds to h/a = —1.91, and (3 = tt/2 + tt/4 = 37r/4, 
which corrrcsponds to h/a = 1.60. By definition, at an 
inflection point dFAo/dh is extremal, which implies that 
the capillary force exerted on the particle by the inter- 
face, defined as 



/» 



dF 



.-id 



Oh 



-f(h), 



(76) 



also has an extremum. However, instead of discussing the 
derivative, which becomes complicated as soon as FAo(h) 
is given implicitly by a set of equations (see Appendix B), 
we resort to a physical picture. 

The total capillary force has two components. The 
first one is the force acting on the contact line at the 
particle, which for a unit length dl of the contact line, 
equals je t dl, where e t is a unit vector tangential to the 
liquid-gas interface and normal to the contact line, point- 
ing in the direction away from the particle (see, e.g., 
Ref. [47j|). In the case of an axisymmetric interface, in- 
tegrating over the contact line leaves only the vertical 
component, equal to — 27T7asin/3sin'0Oj where 27rasin/3 
is the circumference of the three-phase contact line, the 
tangent of the angle —ip(r) equals the slope of the in- 
terface at a distance r from the axis, and, according to 
Young's law, ip = ip(r = asin/3) = j3 — 9 p (see, c.f., 
Appendix B and Fig. ITT]) . The second contribution to 
the total capillary force stems from the internal pressure, 
which for an almost spherical interface is approximately 
equal to the Laplace pressure 2j/Rq (given exactly by 
Eq. (|B15|) ) . resulting in a force of magnitude approxi- 
mately 7r(asin/3) 2 27/i?o directed upwards, which is of 
the order O(a/Ro) relative to the force on the contact line 
and therefore it can be neglected for Ro 3> a. Hence, for 
9 p = tt/2 (as studied in Fig. [5]) the force on the particle 
equals approximately 7T7asin(2/3), with its extremes ap- 
proximately being nja and —irja at j3 = it /A = tt/2— tt/4 
and (3 — 3tt/4 = tt/2 + tt/4, respectively. 

As a consequence of the liquid volume constraint, even 
for large droplets the free energy FAo(h) must be always 
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FIG. 5: The surface free energy FAo{h) (a) for axially sym- 
metric configurations (a = 0) with #o = t/3, v — tt/2, and 
Vi = 79^-a 3 and the total capillary force f(h) = -dF A0 /dh 
(b) as functions of immersion h. For h > and h < one 
can distinguish 7 branches (for reasons of clarity only the 
branches for h > have been tagged in the plot). Branch 1 
corresponds to the globally stable, full analytic solution of the 
Young-Laplace equation. Branch 2 is the free energy under 
the constraint that the droplet keeps the shape of a spher- 
ical cap, which is the equilibrium configuration for h = 0. 
For h ^ this constraint induces 6 V to deviate from its value 
tt/2. Lifting this constraint lowers the free energy towards the 
globally stable branch 1, on which 6 P has its fixed value tt/2. 
For h/a > 0.93 and h/a < —1.04 the particle detaches from 
the interface if one keeps the droplet shape to be a spherical 
cap (3, 4). On branch 3 (4) this detached configuration is 
metastable (globally stable). The extension of branch 1 be- 
yond the horizontal branches 3 and 4 is metastable on branch 
5 up to the filled square, which indicates an inflection point. 
Beyond that point branch 6 is unstable. Branch 7 is an ex- 
tension of branch 2 which is metastable with respect to de- 
tachment, i. e., branches 3 and 4. Branches 6 and 7 merge 
at the point at which the unstable solution along branch 6 
ceases to exist. The equilibrium configurations are given by 
the branches 1 and 4. Lines correspond to analytic results 
whereas circles correspond to numerical finite element calcu- 
lations. The various line codes in (a) and (6) correspond to 
each other. In (b) the straight long-dashed line indicates the 
slope of the analytic solution at h = 0. The force maximum 
fmax/i'jo,) = 3.57 occurs at h/a — —1.82 and the minimum 
/min/(70) = -2-79 at h/a = 1.70. 
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the one with the particle completely immersed in the 
liquid (as long as 9 p = tt/2 the surface energy at the 
particle does not vary with h). Accordingly, the total 
capillary force acting on the particle f(h) is not exactly 
antisymmetric with respect to h = 0. In the studied 
case of a medium-sized droplet the extremal values of 
the force, calculated by a numerical differentiation of 
F A0 (h), equal f ma x/(ja) = 3.57 and /,„;„/ (7a) = -2.79 
for h/a = —1.82 and h/a = 1.70, respectively (see Fig. 
[SJ6)); thus this effect is clearly detectable. From the de- 
viation of the full line from the dashed one in Fig. [5j6) 
one can also infer that the nonlinearity of f(h) is stronger 
for h > 0, i.e., for negative /. The deviation from the 
linear behavior sets in approximately for //(7a) > 2 and 
for //(7a) < -1. 

In Fig. [6] we compare the shapes of the droplet for the 
configurations with (3 = ir/4 (for which f/(ja) = 3.56) 
and /? = 37r/4 (for which //(7a) = —2.77) calculated by 
using the exact solution of the non-linear Young-Laplace 
equation (Appendix B) and from linear perturbation the- 
ory. (For the case of axial symmetry Appendix C pro- 
vides the solution of the linearized Young-Laplace equa- 
tion also for #0 7^ 7 r/2). Apparently, even close to the 
edge of the stability regime the discrepancies between 
the two approaches are so small (actually j3 = ir/4 cor- 
responds already to a metastable configuration, so that 
it belongs to branch 6 for h < in Fig. Eta)), that they 
are practically beyond the resolution of Fig. [SJ we expect 
them to remain small also for the configurations without 
axial symmetry. 



FIG. 6: Droplet shapes for 9 P — tt/2 and fixed 80 = n/3 
(model A), calculated by using the exact solution of the cap- 
illary equation (|B8[) (black solid lines beneath the red dashed 
lines) for /3 = 37r/4 in (a) and for /3 = n/4 in (6), and by using 
the approximate solutions Roevo(0)/a — u(8)/a (red dashed 
lines) given by Eqs. ()C1|) - ()C3|I for pointlike forces of amplitude 
//(7a) = 2.77 and //(7a) = —3.56, respectively (compare 
Fig. [5j6) and Appendix C), with a cut-off at 6 = 0.005 and 
with the remaining set of parameters the same as in Fig. [S] 
We note the logarithmic divergence of the approximate solu- 
tions at 8 — 0. The dotted lines correspond to the particle 
and the spherical cap shapes of the droplet in the reference 
configuration. 



slightly asymmetric with respect to ft, = (or cquiva- 
lently f3 = tt/2 + O(a/Ro)). This can be understood 
by comparing the two most extreme configurations. In 
order to maintain a constant volume the radius of the 
droplet for the configuration with the particle completely 
immersed in the liquid phase (h/a = —1.06 on branch 7, 
see Fig. [5ja)) must be larger than for the configuration 
with the particle completely immersed in the gas phase 
(h/a = 0.94). Because in these two cases the droplet 
shapes are spherical caps, the free energy is higher for 
the configuration with a larger droplet radius, i.e., for 



B. Free energy of configurations without axial 
symmetry 

In order to determine numerically the dependence of 
the free energy on a we fix the angular position a of 
the particle. For such a configuration we minimize the 
free energy functional in Eq. ([3]) with respect to the ra- 
dial displacement h of the particle and thus the shape of 
the droplet by using a finite clement method [46j. This 
nrcthod tracks the evolution of an arbitrarily shaped ini- 
tial body of liquid towards its shape corresponding to the 
minimum of the free energy. In order to obtain suitable 
start configurations with a ^ the surface is pre-evolved 
for a = from a cubic shape and then the particle is 
moved away from the apex in a step-wise fashion by small 
increments 8a. In both model A and model B the con- 
tact line at the particle is taken to be free. In the case 
of model B the contact line at the substrate is pinned 
from the very beginning at a circle corresponding to the 
reference configuration with / = 0; accordingly, the geo- 
metrical center of the reference droplet (see Fig. [1} , from 
which the radial distance of the particle is measured, is 
well defined throughout the evolution. In the case of 
nrodcl A the situation is more complicated because, due 
to the finite size of the particle, the position xcM,ref of 
the center of mass of liquid depends on a. This finite- 
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size effect has to be taken into account for medium-sized 
droplets (see, c.f., the end of Subsec. V.B.2). 
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FIG. 7: (a) Rescaled excess free energy (Eq. (|5J) for a pinned 
contact line on the substrate as a function of the polar an- 
gle a (Fig. [TJ for a contact angle Q Q = 7r/2 of the reference 
configuration at the substrate. The solid black curve cor- 
responds to the analytic expression given by Eqs. (|70|) - (|74|) . 
The dotted line and the dashed line are the contributions 
from the first and the second image (Fig. [2J|, respectively, 
and the dash-dotted line is the volume correction (Eq. (16811 '), 
Green triangles (A), blue circles and red squares correspond 
to {R a /a = 4,//( 7 o) = 2,6 P = tt/2}, {Ro/a = 8,//( 7 a) = 
1,0 P = 2tt/3}, and {Eo/a = 12, f/(ja) = 2,9 P = tt/2}, re- 
spectively. Blue diamonds and inverted triangles (v) core- 
spond to {Ro/a = 8, //(70) = — 1.5, # p = tt/2} and {-Ro/a = 
8, //(7a) = — 2,9 P = 7r/2}, respectively, (fe) The difference 
A9(4>) = 0{(j)) — 0o between the actual contact angle 6 {(f)) and 
6*o = 7r/2 as a function of the azimuthal angle <j) £ [0, it] for 
{Ra/a = 4, //(7a) = 2, # p = 7r/2} and for various angular 
positions a of the particle. The dashed lines correspond to 
the approximate analytic expression (Eq. (J77J) and the solid 
lines are numerical results. The thin horizontal line at A9 = 
corresponds to the reference configuration. 



dependence and we obtain a single master curve, in very 
good agreement with the theoretical expression in Eq. 
([72]) for a pointlike force. For the latter the contributions 
from both images (dashed line and dotted line) and from 
the volume correction Gs,corr (dash-dotted line) are all 
equally important. We emphasize the dc facto indepen- 
dence of AF from the contact angle 9 p at the particle 
(see Fig. 7 and additional data not shown), which jus- 
tifies our perturbation theory which does not take into 
account the specific wetting energy of the particle. The 
only systematic deviation from the analytic expression 
occurs for negative value of /. We have investigated the 
system for various absolute values of negative / (data not 
shown) and the free energy turns out to be always slightly 
overestimated by the analytic theory. In fact we obtain 
two master curves (for fixed Ro/a), which are close to 
each other: one for positive and one for negative /. For 
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FIG. 8: (a) Numerical results for the dependence of the 
rescaled excess free energy (Eq. ®) on the angular position of 
the particle for various values of the contact angle 60 of the ref- 
erence configuration for a pinned contact line on the substrate. 
For all curves f/{-ya) = 2,6 P = tt/2, and Vi = 79 X ^-a z . (b) 
Position and depth of the free energy minimum as functions 
of 9q. The symbols in (6) correspond to those used in (a). In 
(a) the curves end at that value of a at which the particle 
makes contact with the substrate. For pointlike particles one 
has a < 80 . In all plots the lines are guides to the eye. 



1. Pinned contact line on the substrate 

In Fig. UJa) we plot numerical values of jAFb/ f 2 for 
various droplet sizes Rq, both signs and various strengths 
of /, and various contact angles p at the particle. The 
variation of all these parameters does not affect the a- 



/ 7= the actual contact angle 9{4>) at a pinned contact 
line (see Eq. ([341) ) differs from the constant value 80 for 
the reference configuration (/ = 0). Up to first order in 
e one obtains for 9 = tt/2 and with v(fl) — q Gb(Q, ^1) 



cos 



Rq ja 



deGB(Q,Qi)\n=(6=Tr/2,<p)- (77) 



The numerical results for A9(4>) :— 9(<fi) — tt/2 and the 
comparison with Eq. ([77)) are presented in Fig. [7J6). 
The discrepancy is typically of the order of 10%, which 
differs from the almost perfect agreement found in the 
case of the free energy (see Fig. UJa)). If the parti- 
cle is close to the contact line (a = 72°) the discrep- 
ancy reaches 25% which is linked to the large value of 
A6>(</>) w 30° « 0.5 [rad] for cp sa 0, which in turn sig- 
nals that in the vicinity of the contact line and close 
to the particle the small gradient approximation deterio- 
rates and the terms 0(e 2 ) in Eq. (|34p become important 
(for the values of the parameters used in Fig. [7J6) one has 
e = 0.5). Finally, we observe a strong dependence of the 
free energy on do (Fig. [8]) such that upon increasing 9$ 
the free energy minimum moves away from the apex but 
also stays away from the contact line on the substrate 
and its depth increases strongly if 9$ approaches it. 



2. Free contact line on the substrate 

The results for a free contact line, presented in Fig. |H1 
demonstrate that changing the boundary conditions on 
the substrate can change the behavior of the particle 
completely. The global equilibrium position is now at 
the contact line but there is a deep metastablc free energy 
minimum at the drop apex. For forces / of the order of "fa 
the free energy barrier is of the order of 0.05 x 7a 2 . Thus 
for micron-sized particles this is typically much higher 
than the free energy of thermal fluctuations, and there- 
fore the corresponding configuration is expected to be 
experimentally observable. 
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FIG. 9: The rescaled excess free energy (Eq. l[8])) as a func- 
tion of the polar angle a for the contact angle do = tt/2 of a 
free contact line at the substrate, //(7a) = 1 (above the solid 
line) and //(70) = — 1 (below the solid line), and various 
radii Ro. The solid line corresponds to the expression given 
by Eqs. (|65p . (|72p . and (|73p for pointlike forces. For geomet- 
rical reasons the maximal accessible values of a decrease for 
decreasing _Ro- For a free contact line the contact angle at the 
substrate remains at the value 80 for / 7^ 0, i.e., it remains 
the same as for the reference configuration. 

In the case of a free contact line we observe much larger 
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FIG. 10: (a) The rescaled excess free energy (Eq. (|8])) for a 
free contact line at the substrate after taking into account 
the finite size correction SF (Eq. I|79[>) for the numerical data 
(symbols): A' Fa = AFa — SF; the color code and the mean- 
ing of the solid line are the same as in Fig. [9l //(7a) = 1 
and //(70) = — 1 for the points slightly above and below the 
solid line, respectively (visible for a > 60°). (b) The defor- 
mation of a free contact line at the substrate for a = 12° 
(squares), a = 48° (diamonds), and a — 72° (triangles), 
where do = tt/2, Ro/a = 4, //(7a) = 2, and 6 P = tt/2. The 
dashed black line denotes the reference configuration. The 
solid lines are predicted analytically for pointlike forces by 
Eqs. (|49|l . (|65[) , and (|67p . The contact angle is constant along 
these contact lines and equals 80 — tt/2. 



discrepancies with the predictions of the theory for point- 
like forces, revealing a dependence of AFa on / beyond 
the simple scaling ~ f 2 . However, this deviation vanishes 
for increasing radii Rq of the droplet indicating that this 
is a finite-size effect. In order to understand this effect 
we first consider the reference configuration, i.e., the case 
/ = 0. For 6 P = tt/2 the immersed part SV re f of the par- 
ticle (being the intersection of the domain occupied by 
the particle with the spherical cap, representing the ref- 
erence droplet of volume Vi = (27ri?^/3)(l + O(a/i? ) 3 ) in 
the case 9$ = 7r/2) has, independently of a, the volume of 
SVref = (27ra 3 /3)(l + O(a/R )), but the position of this 



18 



cavity in the liquid depends on a. Therefore xcM,ref also 
depends on a and equals (here the position of the particle 
is taken to have a positive x-component, see Fig. Q}: 



XCM,ref(a) 



S v , dV '- 



Vi 

SVref 



1 

v. 



dVx 



i Jsv T ^ 



Rosina~ — — a sin a. (78) 
n \ Rq / 

In contrast, in the case of a pointlike force one has 
xcM,ref — 0. The corresponding difference in the free 
energy SF can be understood as the work done by the 
force fcM = —/sin a (see the main text after Eq. (|43|l ) 
applied to the center of mass in order to counterbalance 
the lateral component of the force / (see Eq. (l43l) ). upon 
displacing the center of mass from the configuration with 
a = to a: 



5F = da' fcM{a')dxcM,ref/da' 
Jo 

t ( a \ r a ' ■ ' i f a ( a \ ■ 2 

= ja — / da sin a cos a = — — - sin a. 
\RoJ Ja 2 \RqJ 

(79) 

This contribution |48[ should be subtracted from the 
numerically calculated free energy in order to facilitate 
the comparison with the analytical result for a pointlikc 
force. It is linear in /, which explains the aforemen- 
tioned deviation from the scaling ~ / 2 . For large drops 
it vanishes ~ (a/Ro) 2 so that for (a/Ro) ^$ 0.1 it can 
practically be neglected (compare Fig. [9]). However, for 
smaller droplets this correction has to be taken into ac- 
count in order to obtain agreement with the perturbation 
theory (see Fig. HPT a)). 



VI. SUMMARY AND CONCLUSIONS 

We have studied the influence of curvature and confine- 
ment on the free energy of a particle of radius a floating at 
the surface of a sessile droplet of radius Rq and exposed 
to external forces / acting in the direction normal to the 
unperturbed droplet surface (see Fig. [T]) . Our results can 
be summarized as follows. 

If the particle is at the drop apex the system exhibits 
axial symmetry which allows one to obtain exact ana- 
lytic solutions for the droplet shape and expressions for 
the surface free energy depending on the immersion of the 
particle into the liquid phase (see Fig. [5]). The compari- 
son of these analytic results with those obtained in Ap- 
pendix C within linear perturbation theory indicates that 
non-linear effects contained in the full Young-Laplace 
equation (see Eq. (|B8[) ) can be neglected even for rel- 
atively large particles (see Fig. [6]). The comparison of 
these full analytic results with those obtained by a nu- 
merical minimization of the free energy (Sec. V) validates 
the high accuracy of the numerical code (Fig. [5]). 



In the cases without axial symmetry the condition of 
balance of forces acting on the droplet in lateral direc- 
tions requires either a fixed lateral position of the cen- 
ter of mass of the droplet (leaving the contact line free: 
model A) or a pinned contact line at the substrate (leav- 
ing the center of mass free: model B). Using a per- 
turbation theory for small deformations of the droplet 
we have derived a free energy functional appropriate for 
both models. In the case of a free contact line at the sub- 
strate, the ensuing changes in the substrate-liquid surface 
energy have been also taken into account (Eq. ©). In 
model A the contact angle at the substrate maintains its 
value 6*o of the unperturbed droplet whereas in model B it 
varies along the contact line and can differ from 0$. The 
pulling of the particle by an external force / has been im- 
plemented by introducing an effective pressure field 7r(51) 
(see Fig. [2] and Eq. (l22l ). which enters the linear Young- 
Laplace equation determining the small deformations of 
the droplet (Eq. (|29l0 . Similarly, the fixing of the cen- 
ter of mass is implemented by an effective pressure field 
kcm(Q) (see Eqs. (|25|) and ([5U|) ). We have shown that in 
the limit of small particles, i.e., for a/Ro — > 0, the free en- 
ergy of the sessile droplet (Eq. (|53[0 . expressed in terms 
of a Green's function satisfying the boundary conditions 
at the substrate corresponding to either a free or a pinned 
contact line (Eqs. (|57]) and ffi§\), does not depend on the 
size of the particle but only on the pulling force / (Eq. 
([60]) '). the contact angle 9$ of the unperturbed droplet at 
the substrate, and on the angular position of the particle 
a (Fig. [TJl. 

If the contact angle Qq at the substrate equals 7r/2 one 
can exploit an analogue of the method of images known 
from electrostatics in order to calculate the surface free 
energy (in excess over the surface free energy of the ref- 
erence configuration of a drop shape given by a spherical 
cap) depending on a. The boundary conditions at the 
substrate are met by introducing an image particle at 
the virtual droplet hemisphere below the substrate sur- 
face. In analogy to electrostatics a pinned contact line 
corresponds to Dirichlet boundary conditions (see Fig. 
I3J and a free contact line with fixed contact angle cor- 
responds to Neumann boundary conditions (see Fig. [3]) . 
A change of the type of boundary conditions leads to a 
change of sign of the capillary charge associated with the 
image particle; accordingly, the excess free energy, which 
is proportional to the interaction energy of the original 
particle with its image, also changes sign. Further anal- 
ysis shows that due to the conditions of force balance 
and volume constraint the Green's function requires ad- 
ditional terms (see Eq. (|M|) ). but they do not change the 
results qualitatively. 

If one treats the particle as a capillary monopole which 
exerts a pointlike force on the fluid interface, the excess 
free energy can be expressed in terms of the Green's func- 
tion G(6) (Eq. (|65])) derived by Morse and Witten [|| 
for a free spherical drop, where 8 is the angle between 
the direction of observation of the interface displacement 
and the direction into which the pointlike external force 
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is applied, both in radial direction from the center of the 
reference configuration, i.e., normal to the surface of the 
unperturbed droplet (see Fig. Q]). For S<1 the Green's 
function exhibits a characteristic logarithmic divergence 
recovering the known corresponding behavior for the case 
of a flat interface; the size of the particle acts as a natural 
cut-off. Compared to this latter case the Green's func- 
tion of the finite-sized drop gives rise to a new interesting 
effect due to the non-monotonic behavior of G(9) for in- 
termediate angles. As a consequence, besides the known 
phenomena of attraction of the particle to the free con- 
tact line and repulsion from the pinned contact line, we 
have found a deep local free energy minimum for the par- 
ticle being at the drop apex for a free contact line (Fig. [9]) 
and a global free energy minimum at an intermediate an- 
gle a (Fig. [T]) for a pinned contact line (Fig. [7]). The 
positions of these free energy minima are de facto inde- 
pendent of the strength / of the external force acting on 
the particle, the droplet radius Rq, and the contact angle 
Op at the particle, while their depths are proportional to 
f 2 . For a particle of radius \\rni and an interface with 
surface tension 7 = 0.057Vm -1 the depth is of the order 
of 10 5 /csT and therefore the corresponding configuration 
is expected to be experimentally observable. A pulling 
force / can be applied to the particle by using optical 
tweezers (see, for example, Ref. [431). This technique 
would also enable one to measure the force tangential to 
the liquid-gas interface corresponding to the derivative 
d(AF a ) J d(Roa) of the free energy (a = A, B) with re- 
spect to the lateral displacement of the particle along the 
interface. 

Moreover, our numerical results in Fig. [5] show that for 
model B the depth of the free energy minimum increases 
rapidly upon increasing the contact angle 9q at the sub- 
strate towards 180°. In the case of model A, the shape 
of the contact line at the substrate as obtained from nu- 
merical calculations agrees very well with the predictions 
of the analytic expressions given in Eqs. (|67|) and (|65j) 
(see Fig.HUTb)). whereas the free energy (Fig. [9]), due to 
the condition of fixed center of mass, is more sensitive 
to the finite size a of the particle and converges to the 
expression for a pointlike particle for a/Ro < 0.1. On 
the other hand, in the case of a pinned contact line the 
finite-size effects are negligible (Fig. Eta)). In this case 
the contact angle differs from its reference value 9q and 
varies along the contact line (Fig. EJo)), in agreement 
with the analytic approximation (Eqs. (|34| and (f77| for 
On = tt/2). 

The overall very good agreement with the pointlike 
force approximation, even for relatively large particles 
such as with a/Ro = 0.25, indicates that higher order 
corrections in a/Ro are very small. Indeed, in the case 
of a spherical particle the capillary dipole must vanish 
due to the vanishing torque on the particle and the first 
sublcading term stems from the interaction energy of a 
monopolc with an induced quadrupolc yielding a correc- 
tion of the order (a/i?o) 4 - 
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Appendix A: Free energy of a sessile droplet 

In this appendix we derive the expression for the free 
energy of a sessile droplet in mechanical equilibrium and 
subjected to a pressure field 7r(J7) (Eq. (|53|) ). In terms 
of perturbation theory, the free energy both for model 
A and model B is given by the functional in Eq. (|2"5)) 
evaluated for v obeying the Young-Laplace equation (|2"9")l 
with -kcm = for model B and with the condition in Eq. 
for model A. This yields 



F 



— I 



7 Jn Q 



-V a (wV a w) 



2^(Va^ 



2v) 



■(ir(Q) + n)v] - — cos6>o / #(t>|e ) 

— / d<f)v\e (sin d d, 
^7 Jo 



'gv\ So - cos Q v\e o ), (Al) 



where the second equality follows from applying Gaufi' 
theorem, Eq. (|29|) . Eq. (|36|) . and the constant volume 
constraint (Eq. (|35p ). The second term in the last ex- 
pression vanishes for boundary conditions corresponding 
to either a free or a pinned contact line at the substrate 
(see Eqs. (|32|) and (|33jl). This leads to the free energy in 
the form in Eq. ([55)1 . 



Appendix B: Exact results for axially symmetric 

configurations with a free contact line at the 

substrate 



In this appendix we study the shape and the free en- 
ergy of the droplet with the adsorbed particle at its apex, 
resulting in an axially symmetric configuration. In the 
following calculations we take a as the unit of length and 
7a 2 as the unit of energy. For simplicity we assume that 
the interface has no overhangs (9q < tt/2) so that the 
shape of the droplet can be described by z = z(r). We 
also assume a free contact line and a given contact angle 
9q at the substrate (see Fig. [TTj) . A similar analysis can 
be performed also for a pinned contact line, in which case 
the surface energy associated with the substrate would be 
constant. The problem with a free contact line is actually 
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FIG. 11: Cross section of an axially symmetric configuration 
of the droplet and the particle; (a) reference configuration for 
/ = 0, which leads to a liquid-gas interface with the shape 
of a spherical cap; (6) sketch of a configuration for f3 ^ /3o 
corresponding to / =^ 0. This sketch corresponds to model A 
for which 0q is kept fixed and r m 7^ _Ro sin #0 . The dashed lines 
denote continuations of the analytical profile corresponding to 
an unduloid or a nodoid beyond the physical regime bounded 
by r = sin/3 and r = r m and terminating at maximal and 
minimal radii tq and n, respectively, with z(ro) — z(ri) = 
-co. Note that even for h > 1 this configuration can be stable 
(compare Fig. E£a)). 



analytically more involved, because in this case the sur- 
face free energy associated with the substrate can vary. 
Therefore it is particularly interesting and instructive to 
incorporate the contribution from the liquid-substrate in- 
terface into the total free energy. Our goal is to calculate 
the free energy i^o, introduced implicitly as the inverse 
Legendre transform of Fao in Eq. (|75|) (see also Eqs. © 
and ©) with a = A corresponding to a free contact line, 



as a function of the vertical position h of the particle and 
the independent system parameters #0, p , Rq (the latter 
expressed in units of a). To this end we first consider 
the free energy (compare the analogous expression in Eq. 
(|4|, noting that here there is no force /cm needed): 

FXo(h,6 ,e p ,R ,X) - min T[{z(r)}], (Bl) 
{>M} 

with the free energy functional T as 



f[{z(r)}; h, 6»o, p ,V h A] = S ig - Si g>re f 

— (Soi — SoLref) COS 6>o — {S p i — S p i^ re f) COs8 } 



■X(V-Vi), 

(B2) 



where the volume Vi of liquid is a function of the contact 
angles 9 and 8 P , as well as of the droplet radius R 
in the reference configuration, i.e., Vi = Vi(0o,9 p ,Ro). 
The Lagrange multiplier A can be determined from the 
condition V = Vi, which renders A = \(h,9o,8 p ,Ro) and 
which upon insertion into F\ (h, 6q,0 p ,Rq^ A) yields the 
desired free energy FAo(h, 9 , P , Ro)- 

It is convenient to express the contact areas Soz, S p i, 
and Si g in terms of two auxiliary variables: the angle /3 
describing the position of the contact line at the particle 
and the radius r m of the circle formed by the contact 
line at the substrate (see Fig. [PTT fc)). In terms of these 
variables one has 



Soi=Kr rn , (B3) 

S p i = 2vr(l + cos^) = 4ttcos 2 (/3/2), (B4) 

Sig = 27T 



drr\/l + z' z 



(B5) 



sin j3 



with z = dz{r)/dr. The volume of liquid V can be ex- 
pressed as 

Ajr 

V = Vi g + 7r(sin 2 (3)z(r = sin/3) - — / (tt - /?) (B6) 
with fo(x) given by Eq. (|2"5)) and 



Via 



2tt 



dr rz(r), 



(B7) 



sin /3 



where Vi g is the volume enclosed between the substrate 
and the liquid-gas interface, the second term in Eq. (|B6[) 
is the volume between the substrate and the circular disc 
determined by the contact line at the particle, and the 
third term is the volume of that part of the particle which 
is immersed in the liquid. 

The equilibrium profile z(r), which minimizes T obeys 
the Young-Laplace equation in cylindrical coordinates 
(see, e.g., Ref. (H): 



rz 



1 d 

r dr y/l + z 2 



-X. 



(B8) 



The boundary conditions are determined by Young's law 
at the particle and at the substrate which can be shown 
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to follow from the condition of vanishing of the variation 
of T at the boundaries: 



V>(r = sin 0)= f3 -0 p , 
il>(r = r m ) = 6 , 



(B9) 

(BIO) 



where ip = ip(r) is the angle between the r-axis and the 
tangent of the profile z(r) (see Fig. HH&)) defined by 



sin0(r) 



The first integral of Eq. (|B8[) reads 



■ i , \ Xr c 

■ am ip{r) = — - + -, 
2 r 



(Bll) 



Vl + i 



rl 



(B12) 



where c is an integration constant. When evaluated at 
the boundaries and using Eqs. (|B9j) and (jBlOj) this leads 
to the following set of equations: 



sin(/3 



Asin/3 



sin/3' 



Xr„ 



sin 0o = — " 



2 r m 
Solving with respect to A and c one obtains 



(B13) 
(B14) 



2[r m sin6>p - sin/3sin(/3 - 6 p )} 

. o : =K r m,P), (B15) 



r m - sin P 

r 2 ^ sin /3 sin(/3 — p ) — r m sin Op sin /3 



c(r m ,/3). 
(B16) 



The equilibrium profile z(r) is obtained by solving Eq. 
(|B12[) for z and by subsequently integrating: 



z(r) 



dr- 



Xr 2 - 2c 



Aisin/3 V(r 2 -r 2 )(r 2 -r 2 ) 



2c 



-r [£(0, g) - E(&, g)] + — [F(0, q) - F(0 2 , g)] • 

(B17) 



The integral in Eq. (|B5[) can be evaluated as 

47T 



^S 



<ir- 



|A| y sin /3 \/F^)(^) 



= 4 7 r^LE(g)- J E(0i,g)- J E(0 2 ,g) 

+ ^n7V / ( r o-^ 2 /3)( S m 2 /3-^)). (B18) 

The expression in Eq. (|B7|) can be evaluated by integrat- 



ing by parts, which leads to 

Ar 4 - 2cr 2 
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dr- 
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(r 2 , -sin 2 £)(sin 2 /3-r 2 ) 



- r -Y\J{rl-r 2 m ){r 2 m -r 2 )\ - 7r(sin 2 0)«(r = sin/?), 

(B19) 

where -F (0, fc) and i?(</>, fc) are the incomplete elliptic in- 
tegrals of the first and second kind, respectively, whereas 
K(k) and E(k) are the corresponding complete elliptic 
integrals (see Ref. j51(). Moreover, one has 
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and 
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(B20) 
(B21) 

(B22) 
(B23) 



The shape of the droplet is always a section of a nodoid 
or an unduloid [50( characterized by a maximal and a 
minimal radius rp and n (see Fig. I11T 6)). respectively, 
which are two distinct solutions of the equation z(r) — 
-co and are given by 



1 



ro = j ( VT 1 ) , 



ri = - ( VI 



(B24) 
(B25) 



So far the interface profile (Eq. (|B17|0 and the free 
energy (Eqs. (|Bl|l - ((B4)l and (|B18|l ) are determined in 
terms of the independent variables r m and /3 (through 
Eqs. (|BT5|) , (JBl"6| and (|B20)l - (|B25)l ). One can replace 
them by the physically more directly accessible variables 
h and V;. To this end we note that the vertical displace- 
ment h of the particle is determined by (see Fig.QTJfo)) 



z + h — z(r = sin/3) — cos/3, 



(B26) 



where z — Zq(6o,6 p , Rq) is the height of the particle 
center above the substrate in the reference configuration 
with the droplet shape given by the cap of a sphere (see, 
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c.f., Eq. (|B32I) and Fig.[l2]). This renders the set of equa- 
tions 

h = h(r m ,0,Vi), (B27) 

V l = V(r mi /3), (B28) 

where h(r m , (3, V{) is given by Eq. (|B26|) (with the depen- 
dence on r m entering via z(r) and the dependence on Vi 
entering via z (R ) and Eq. (|B29[) 1 and V(r m ,{3) is de- 
fined (suppressing the explicit dependence on O an d 9 P ) 



via E q s. (|B6| and (|B19| toget her w ith (|B15|) . (|B16| and 
(|B20|) - (|B25l) . Equations (|B27|) and ([B28| provide implic- 
itly the maximal radius r m and the angle (3 as functions 
of Vi and /i. 




FIG. 12: Geometry of the reference configuration. The angle 
Po can be found by comparing two expressions for the radius 
of the contact line at the particle: asin/3o = RosinS, where 
8 can be found to be equal to /3o — 8 P (which then yields 
Eq. ()B30[) ). For calculating the volume Vi of the liquid one 
must subtract from the volume of the spherical cap of radius 
Ro the volume of the immersed part of the particle (beneath 
the contact line CL) and in addition the small spherical cap 
indicated as a hatched region (see Eq. ([B29J). 

Finally, with keeping in mind that here all length scales 
are measured in units of a, the constant liquid volume Vi 
can be expressed in terms of the materials parameters 9$ 
and 9 p as well as the drop size Rq . As one can infer from 
the geometrical features shown in Fig. [12l Vi is given by 



Vi(Ro,9 o ,0 p ) = y([/o(0o) - io(/3 - 9 p )]Rl 

-/o(tt-A))), (B29) 



with 



so that 



/3 (R ,9p) = arcsin (R sin 9 P )/(R„ - 2R cos 9 P + 1) 



1/2 



(B31) 



and with the geometric factor fo(x) given by Eq. (|28|) . 
Additionally, since the presence of the particle does not 
change the contact angle 00, the liquid volume Vi can be 
expressed in terms of the size -Ro = [3V//(47r/ (6'o))] 1 / 3 
of the droplet without the particle. The vertical position 
zq of the particle in the reference configuration can be 
expressed as 

z (Ro, #o, 9 P ) = Ro cos(/3 - 9 P ) - R cos 9 - cos (3 . 

(B32) 
The sets of equations (|B27j) - (JB32J) (which must be 
solved numerically) implicitly yield j3 = /3(/i, Rq) and 
r m = r m (h,Ro), which can be substituted into the ex- 
pressions for the surface areas in Eqs. (|B3[) . (|B4|) . and 
(|B18p . Finally, the free energy functional T given as a 
combination of these areas (Eq. (|B2[) ) yields the free en- 
ergy FAo(h, 9q,9 p ,Ro). The results of these calculations 
for 9 = tt/3, 9 p = tt/2, and V t = 79 X (4tt/3) are pre- 
sented in Figs. [5] and |B1 



Appendix C: Droplet shapes from perturbation 

theory for configurations with axial symmetry and 

arbitrary 8q 

A modification of the method of images can be used for 
arbitrary 9q if the point force / = qe^/Ro is placed at the 
apex of the sessile droplet. The virtual reference droplet 
(for / = 0) is constructed as a smooth continuation of the 
actual droplet such that it completes the full sphere and 
the image of amplitude /' = / is placed at the "south 
pole" 9' = it. In this configuration forces acting on the 
union of the real and the virtual parts of the droplet 
are balanced. We exploit the freedom to add a constant 
and a term proportional to cos 9, such that the axially 
symmetric solution 

v (9) = q[G{9) + G(tt -9)+H cos 9 + I„] (CI) 

conserves the volume and fulfills the given boundary con- 
ditions. Moreover, this solution fulfills Eq. ((29)) with 
7r(f2) = q8(fl) and ttcm(Q) = upon choosing Iq = 
-At/ 2 - 1/(4"") ■ The value of /i, given by Eqs. (gSJ) and 
(|4"?Jf for a = 0, depends on the boundary condition at 
the substrate. In order to be able to make a comparison 
with the nonlinear theory in Appendix B here we take 
the contact line at the substrate to be free. (For the case 
of a pinned contact line see Appendix B in Ref. |23|.) 
Using the expression in Eq. (j4"5j) for a = yields 
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Inserting Eqs. (|C1[) and (|C2[) into the volume constraint 
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Due to I Q (0o = tt/2) = = H (d = ir/2), for O = tt/2 
the solution in Eq. (jCll) reduces to 

w o (0; 6 = tt/2) = g[G(0) + G(7r-0)]sgGU (Q,Oi = 0), 

(C4) 
which reproduces the result of the perturbation theory 
for model A (Eq. (jrJT)) ) for a = 0. In order to express the 
deformation it in the dimension of length the solution in 
Eq. (|Clj) has to be multiplied by Roe = I/I/7. The value 



of the external force / should be taken as minus the capil- 
lary force / (Eq. (|7d| ) acting at the particle. Indeed, due 
to J- = J- + fh, the equilibrium condition dF /dh\ eq = 
is equivalent to / = OF /dh\ eq = dFAo/dh = —f. The 
capillary force f(h) is calculated by numerical differen- 
tiation of the free energy Fa (Ii). Particularly, by us- 
ing the relation h(/3) (Eq. (|B26[) ). one obtains the val- 
ues //(7a) = -3.56 and //(7a) = 2.77, for j3 = tt/4 
and j3 = 37r/4, respectively, for which the curves z(r), 
given implicitely by z(6) = [Rq + u(9)]cos9 — z su bstr 
and r{6) = [Rq + u{6)} sin 6 with u{6) = Roev (6) and 
Zsubstr = Ro cos$o, have been plotted in Fig. [5] 



[1] P. Pieranski, Phys. Rev. Lett. 45, 569 (1980). [22 

[2] J. M. Kosterlitz and D. J. Thouless, J. Phys. C: Solid 

State Phys. 6, 1181 (1973). [23 

[3] K. Zahn and G. Maret, Phys. Rev. Lett. 85, 3656 (2000). 
[4] R. Bubeck, C. Bechinger, S. Neser, and P. Leiderer, Phys. [24 

Rev. Lett. 82, 3364 (1999). 
[5] N. Bowden, A. Terfort, J. Carbeck, and G. M. White- [25 

sides, Science 276, 233 (1997). 
[6] L. E. Helseth, R. M. Muruganathan, Y. Zhang, and T. M. [26 
Fischer, Langmuir 21, 7271 (2005). [27 

[7] J. Aizenberg, P. V. Braun, and P. Wiltzius, Phys. Rev. 

Lett. 84, 2997 (2000). [28 

[8] J. C. Loudet and B. Pouligny, EPL 85, 28003 (2009). 
[9] S. U. Pickering, Journal of Chemical Society 91, 2001 [29 
(1907). 

[10] A. D. Dinsmore, M. F. Hsu, M. G. Nikolaides, M. Mar- [30 
quez, A. R. Bausch, and D. A. Weitz, Science 298, 1006 [31 
(2002). 
[11] M. Chavez-Paez, P. Gonzalez-Mozuelos, M. Medina- [32 
Noyola, and J. M. Mendez-Alcaraz, J. Chem. Phys. 119, 
7461 (2003). [33 

[12] P. X. Viveros-Mendez, J. M. Mendez-Alcaraz, and 

P. Gonzalez-Mozuelosa, J. Chem. Phys. 128, 014701 [34 

(2008). [35 

[13] A. R. Bausch, M. J. Bowick, A. Cacciuto, A. D. Dins- [36 
more, M. F. Hsu, D. R. Nelson, M. G. Nikolaides, 
A. Travesset, and D. A. Weitz, Science 299, 1716 (2003). [37 

[14] J. Ruiz-Garcia, R. Gamez-Corrales, and B. I. Ivlev, Phys. 

Rev. E 58, 660 (1998). [38 

[15] F. Ghezzi and J. C. Earnshaw, J. Phys.: Condens. Matter 

9, L517 (1997). [39 

[16] F. Ghezzi, J. C. Earnshaw, M. Finnis, and M. McCluney, 

J. Colloid Interface Sci. 238, 433 (2001). [40 

[17] R. P. Sear, S. W. Chung, G. Markovich, W. M. Gelbart, 

and J. R. Heath, Phys. Rev. E 59, R6255 (1999). [41 

[18] M. G. Nikolaides, A. R. Bausch, M. F. Hsu, A. D. Dins- 
more, M. P. Brenner, D. A. Weitz, and C. Gay, Nature 
420, 299 (2002). [42 

[19] M. M. Nicolson, Proceedings of the Cambridge Philo- [43 

sophical Society 45, 288 (1949). 
[20] P. A. Kralchevsky, V. N. Paunov, I. B. Ivanov, and K. Na- [44 
gayama, J. Colloid Interface Sci. 151, 79 (1992). [45 

[21] D. Stamou, C. Duschl, and D. Johannsmann, Phys. Rev. 

E 62, 5263 (2000). [46] 



J. C. Loudet, A. M. Alsayed, J. Zhang, and A. G. Yodh, 
Phys. Rev. Lett. 94, 018301 (2005). 

M. Oettel, A. Dommguez, and S. Dietrich, Phys. Rev. E 
71, 051401 (2005). 

A. Dommguez, M. Oettel, and S. Dietrich, J. Chem. 
Phys. 127, 204706 (2007). 

H. Lehle, E. Noruzifar, and M. Oettel, Eur. Phys. J. E 
26, 151 (2008). 

M. Megens and J. Aizenberg, Nature 424, 1014 (2003). 
L. Foret and A. Wiirger, Phys. Rev. Lett. 92, 058302 
(2004). 

A. Dommguez, M. Oettel, and S. Dietrich, J. Phys.: Con- 
dens. Matter 17, S3387 (2005). 

K. Danov and P. Kralchevsky, Adv. Colloid Interface Sci. 
154, 91 (2010). 
M. Oettel and S. Dietrich, Langmuir 24, 1425 (2008). 

A. Dommguez, M. Oettel, and S. Dietrich, J. Chem. 
Phys. 128, 114904 (2008). 

H. Diamant, Journal of Physical Society of Japan 78, 
041002 (2009). 

B. X. Cui, H. Diamant, and B. H. Lin, Phys. Rev. Lett. 
89, 188302 (2002). 

A. Wiirger, EPL 75, 978 (2006). 

A. Wiirger, Phys. Rev. E 74, 041402 (2006). 

A. Dommguez, M. Oettel, and S. Dietrich, EPL 77, 

68002 (2007). 

P. A. Kralchevsky, V. N. Paunov, and K. Nagayama, J. 

Fluid Mech. 299, 105 (1995). 

P. A. Kralchevsky and K. Nagayama, Particles at Fluid 

Interfaces (Elsevier, Amsterdam, 2001). 

A. Sangani, C. Lu, K. Su, and J. Schwarz, Phys. Rev. E 

80, 011603 (2009). 

L. Schimmele, M. Napiorkowski, and S. Dietrich, J. 

Chem. Phys. 127, (2007). 

A. Dommguez, in Structure and functional properties of 

colloidal systems, edited by R. Hidalgo-Alvarez (CRC 

Press, Boca Raton, 2010), pp. 31-59. 

A. Dommguez, private communication. 

J. D. Jackson, Classical Electrodynamics (Wiley, New 

York, 1975), 2nd ed. 

J. Guzowski, PhD thesis, unpublished. 

D. C. Morse and T. A. Witten, Europhys. Lett. 22, 549 

(1993). 

K. Brakke, Experimental Mathematics 1, 141165 (1992). 



24 



[47] L. A. Segel, Mathematics applied to Continuum Mechan- 
ics (Dover, New York, 1987). 

[48] In calculating 8F we have ignored the correction Sx, 
which also depends on a (see Eq. (JED, but gives a con- 
tribution of the order {f 2 /-y) x O((a/R ) 3 ). 

[49] V. Blickle, J. Mehl, and C. Bechinger, Phys. Rev. E 79, 
060104 (2009). 



[50] D. Langbein, Capillary Surfaces (Springer, Berlin, 2002), 

2nd ed. 
[51] A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev, 

Integrals and Series, Vol. 1 (Gordon and Breach, New 

York, 1986), 2nd ed. 



